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Abstract. A system of N non-canonical dynamically free 3D harmonic oscil- 
lators is studied. The position and the momentum operators (PM-operators) 
of the system do not satisfy the canonical commutation relations (CCRs). In- 
stead they obey the weaker postulates for the oscillator to be a Wigner quan- 
tum system. In particular the PM-operators fulfil the main postulate, which 
is due to Wigner: they satisfy the equations of motion (the Hamiltonian's 
equations) and the Heisenberg equations. One of the relevant features is that 
the coordinate (the momentum) operators do not commute, but instead their 
squares do commute. As a result the space structure of the basis states corre- 
sponds to pictures when each oscillating particle is measured to occupy with 
equal probability only finite number of points, typically the eight vertices of 
a parallelepiped. The state spaces are finite- dimensional, the spectrum of the 
energy is finite with equally spaced energy levels. An essentially new fea- 
ture is that the angular momenta of all particles are aligned. Therefore there 
exists a strong interaction or correlation between the particles, which is not 
of dynamical, but of statistical origin. Another relevant feature is that the 
standard deviations of, say, the A;th coordinate and the momenta of ath is 
ARak^Pak < ph/\N — 3| { N ^ 3, p— fixed positive integer), namely instead 
of uncertainty relations one has "certainty" relations. The underlying Lie su- 
peralgebraic structure of the oscillator is also relevant and will be explained in 
the context. 
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1. Introduction 

The title of the present paper comes to indicate from the very beginning that the geometry 
of the quantum system under consideration is noncommutative. To be more exphcit, we 
study A?"— particle three-dimensional harmonic oscillators with a Hamiltonian 

a=l 

such that the position operators -Rai, Ra2^ RaZ do not commute with each other and the 
momentum operators -PqiI: Pa2i Pa3 do not commute, too: 

[Rar, R^j] 7^ 0, 7^0, = 1,2,3. a, (3 ^ 1, 2, N. (1.2) 

For this reason the oscillator system is more involved to study. On the other hand however 
it is simpler, because the squares of all position and momentum operators (PM-operators) 
do commute with each other, 

[Rl„Pl]=0, [Rl„R%] = 0, [P2.,p|.]=0, = 1,2,3. «,/?=l,2,...,Ar, (1.3) 

which is not the case for the canonical oscillator. 

The motivation for such "commutation" relations will be clear soon. Here we remark 
only that the above relations are not postulated. We derive them. 

The other part of the title, ferromagnetic-like oscillators, is to stress on the very strong 
statistical interaction between the angular momentums of the oscillating particles. Despite 
of the lack of any dynamical interaction term in the Hamiltonian, all angular momentums 
of the particles are aligned, they point into one and the same direction, similar to the spins 
in ferromagnets. In the present case however the particles are spinless and they carry no 
electric charge. 

Three other appropriate " candidates" for a place in the title were: 
1. "A Lie super algebraic approach to quantum statistics", coming to indicate that 
with each infinite class of basic Lie superalgebras A, B, C or V one can associate quantum 
statistics, namely particular for this class "commutation" relations between the position 
and the momentum operators (PM-operators). In this terminology the Bose and the Fermi 
statistics are B statistics, whereas the statistics of the present model is, as we shall see, A 
statistics and more precisely sl{3\N) statistics. 
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2. "Discrete quantum systems" pointing out that (as a rule) the space structure of 
most of the basis states corresponds to a picture when the oscillating particle is measured 
to occupy with equal probability only finite number of points, typically the eight vertices of 
a parallelepiped (see figure 1, p. 35). This is another strong correlation property, because 
in most cases these points are the same for all particles. 

3. "Finite-level quantum systems" , a small subtitle indicating that the results of the 
present paper can be of interest also in the context of quantum computing. 

There are several other properties, which differ from those of canonical oscillators 
with Hamiltonian (1.1). Some of them are evident as for instance that there exist neither 
coordinate nor momentum representation. Other properties are not so evident. One of 
them, which will be derived in the present paper, is the analogue (or, rather, an "anti"- 
analogue) of the uncertainty relations. It reads: 

A^«feAP«fe < P^^, iV^3, (1.4) 

where p is a fixed positive integer, labelling the state space under consideration. The 
above inequality holds simultaneously for any particle a and any coordinate k. Note that 
contrary to the canonical case here the left hand side is smaller than the right hand side. 

Passing to a more systematic exposition, we recall the definition of a Wigner quantum 
system. It is based on the following six postulates [1] (in the Heisenberg picture for 
definiteness; the very name WQS was introduced in [2]): 

(PI) The state space W is a Hilbert space. To every state of the system there 
corresponds a normed to 1 vector from W . 

(P2) To every physical observable L there corresponds a Hermitian (self-adjoint 
and hence linear) operator L in W. 

(P3) Given a physical observable L, the measurement outcome values it may assume, 
are just the eigenvalues of the operator L. 

(P4) The expectation value of L in a state ij) is given by (L)* = (V', Lij)). 

(P5) Ri, . . . , Riv and Pi, . . . , Pjv are solutions of the equations of motion 
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(the Hamilton's equations), which for the Hamiltonian (1.1) read: 




R, 



■a 




for a = 1, . . . , A?". 



(1.5) 



(P6) Ri, . . . , Rjv and Pi, . . . , Pat are solutions of the Heisenberg equations 




[iy,Pa], R, 



■a — 



n 



[iy,Ra], for a = l,...,Ar. 



(1.6) 



A WQS with harmonic oscillator Hamiltonian (and in particular the Hamiltonian 
(1.1)) is said also to be (A'"-particle 3D) Wigner quantum oscillator (WQO) or simply 
Wigner oscillator. 

The only difference between the above postulates and the postulates of conventional 
quantum mechanics is that in the latter the postulate (P6) is replaced with the 

(P6) Ri, . . . , Rjv and Pi, . . . , Pat, satisfy the canonical commutation relations (CCRs) 
(see [3], we wrote them as given in [4]): 



[Rai,P/3j] = i-lT'Sai,l3j, [Rai, R/Sj] = [Pai, P/Sj] = 0, i,j = 1,2,3. a, P = 1,2, N. 



Postulates (PI) - (P5) remain the same. 

A necessary step on the way to establish whether a given Hamiltonian admits an alter- 
native statistics, i.e. admits alternative "commutation relations" between PM-operators, 
is to solve the so called Wigner 's problem, namely to find out whether the postulates 
(P5) and (P6) admit common noncanonical solutions. In our case this means to solve the 
compatibility conditions 



Our main task is to find (noncanonical) solutions of the above equations. We shall 
come back to this problem in the next section. At this place we postpone in order to say 
a few words as a justification for the definition of a WQS given above. In doing so one has 
to answer to at least three questions: 

(a) Why to replace the postulate (P6) with (P6)? 

(b) If so, does one obtain in this way new, different from (1.7), relations among the 
position and momentum operators? 



(1.7) 



[H, Pa] = ihniaUj'^'Ra, [H , Rq] 




for a = 1, . . . , N. 



(1.8) 
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(c) Do the new relations lead to new and interesting predictions? 

The first two questions were actually raised and answered by Wigner in a two page 
publication [5] in 1950. First of all Wigner gave a positive answer to question (b). To this 
end he considered an example of a one-dimensional harmonic oscillator with a Hamiltonian 
{m = Lij = h = l)H = |(-P^ + R?). Abandoning the canonical commutation relations 
(CCRs) [P, R] — —i, Wigner was searching for all operators R and P, such that the 
"classical" equations of motion P = —R, R = P were identical with the Heisenbcrg 
equations P = —i[P, H], R = —i[R, H]. The result: Wigner found infinitely many solutions 
labelled by one positive integer p = 1,2, .... Only the p = 1 solution coincided with the 
canonical R and P. In other words Wigner has shown that the CCRs can be viewed 
as sufficient, but not necessary conditions for the Hamiltonian equations and Heisenbcrg 
equations to hold simultaneously. The canonical postulates (PI) - (P5) and (-P6) imply 
that any conventional quantum system is a Wigner quantum system, but they do not 
exhaust all WQSs. 

Next Wigner answered also the question (a) noting that the Heisenbcrg equations 
(1.6) and the Hamilton's equations (1.5) have a more immediate physical significance than 
the CCRs. Therefore it is logically justified to postulate from the very beginning these 
equations instead of the CCRs (1.7). 

Turning to question (c) we list some of the characteristics of WQSs studied so far. 
WQSs from the class A [6] basic Lie superalgebras [1, 7-10]: 

(i) The state space is finite-dimensional. 

(ii) The spectrum of the energy is equally spaced but finite. 

(iii) The geometry is noncommutative. 

(iv) The spectrum of the position operators is finite. 

Very different are the properties of the WQSs related to the LSs from the class B. In 
particular for the WQS related to the LS osp{3/2) from this class [11,12]: 

(v) The orbital momentum of two spinless particles curling around each other 
can be 1/2. This would mean that the spin has a classical analogue. 

(vi) The state space is infinite-dimensional. 

Various aspects of Wigner's idea were studied by several authors from different points of 
view. Among the earlier papers we mention [13 - 18], but the subject is of interest also now 
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[19 - 27]. Here are in short some of the results. In [14] Schweber extended the conclusions 
of Wigner to QFT showing (on an example of a scalar field) that the field's commutation 
relations are also not defined uniquely. Okubo [17] related the different solutions of the 
Wigner's problem (different quantization) to the circumstance that different Lagrangian 
may lead to one and the same equation of motion. In [19] the Wigner's problem was solved 
for a magnetic dipole precessing in a magnetic field, thus demonstrating that the equations 
of motion can be compatible with the Heisenberg equations not only for potentials of 
oscillating type. It is particularly interesting also that the class of noncanonical solutions 
determined in [19] includes the deformed CCRs [28], thus indicating that the quantum 
deformations can be viewed also as generalizations of quantum statistics. 

A strong "push" for studying further alternative commutation relations came also from 
the predictions of string theory that the geometry of the space becomes noncommutative at 
very small distances (see [29] for a survey and the references therein). To similar conclusions 
lead also various deformed models (most of them in the sense of quantum groups (see [30] 
for a review and the references therein). However, the idea itself was already suggested by 
Heisenberg in the late 1930's (as explained in [31]) and perhaps the first example of this 
kind was given by Snyder [32] . 

The paper is organized as follows. In the next sections we recall shortly where the idea 
for application of the Lie superalgebras in quantum statistics comes from. This section 
contains no new results. 

In section 3 we outline the mathematical structure of the sl{3\N) WQO. We identify 
the position and the momentum operators Ri, . . . , Rjv and Pi, . . . , Pjv as odd operators 
in such a way that the linear span of these operators and their anticommutators close the 
Lie supcralgcbra s/(3|A^). Already in this chapter it becomes evident that the angular 
momentum operators poses unusual properties. Up to multiplicative constants they are 
the same for all particles and coincide with the angular momentum operators for the entire 
system. Mathematically this means that the projections of the angular momentum of the 
different particles are described by operators, which are equal in the sense of operators (up 
to multiplicative constants). Physically it corresponds to a picture with angular momentum 
alignment of all particles. 

In section 4 the state spaces V{N,p), p = 1,2, 3, ... of the WQO are introduced. These 
are finite-dimensional subspaces of the infinite-dimensional Fock space W{Z\N), generated 
by three pairs of Bose creation and annihilation operators (CAOs) and N pairs of Fermi 
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CAOs. Each V{N,p) carries an irreducible representation of sl{3\N). The state spaces 
corresponding to different p carry inequivalent representation of the LS sl{3\N). A some- 
what unusual feature of this construction is that the Bose operators are considered as odd 
generators, the Fermi operators are even elements and the Bose operators anticommute 
with the Fermi operators [33]. This unconventional grading of the Bose and Fermi op- 
erators is not accidental. In this way W^(3|iV) carries an infinite-dimensional irreducible 
representation of the orthosymplectic LS osp{3\N). 

In section 5 the energy spectrum of the entire system and of each of the oscillating 
particles is derived. The energy spectrum of the system is equidistant but finite. The 
Hamiltonian (1.1) has min{N,p) + 1 equally spaced energy levels with a gap between 
neighboring levels uh. The multiplicity of each such level is computed. The ground state 
is nondegenerate only in the case when N = p. In the case N = 1, namely for one 3D 
WQO the energy levels are only two and if in addition p = 1 they coincide with the first 
two energy levels of a canonical 3D oscillator. 

The energy of each individual oscillator is an integral of motion. The energy levels 
are again equidistant, but this time the gap between neighboring levels is a fraction of coh 
and more precisely it is u)?i/\N — 3|. One of the unexpected features of the single particle 
energy is that in certain cases its ground energy can be zero (together with coordinates, 
momenta and angular momentum). 

Section 6 is the biggest one. Here the space structure of each basis vector from the 
Fock space is analyzed. It is shown that certain states correspond to a picture when a 
particle can be measured to be with equal probability on every point of a sphere. There 
are also states with a particle distribution along two circles (figure 5, p. 69), but the typical 
picture is that the particle is measured to occupy with equal probability the eight vertexes 
of a parallelepiped (see figure 1, p. 34). For any Fock state the standard deviation of the 
particles along any direction is computed too. 

In Section 7 the so(3) structure of the state space V{N,p) is clarified. Each such 
space is decomposed into irreducible so(3) modules. As already mentioned, up to overall 
constants the projections of angular momentums are the same for all particles, which 
results in the angular momentum alignment. In this section the parity operator of each 
particle is introduced and the property that the nests of the basis states are occupied with 
equal probability is proved. 

The last section 8 contains some concluding remarks. 
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Some abbreviations and notation: 



[x, y]=xy- yx, {x, y} = xy + yx. 

{(p,ip) - the scalar product between the states (p and ip; 

Z+ - all nonnegative integers, 

N - all positive integers, 

WQS(s) - Wigner quantum system (s), 

WQO(s) - Wigner quantum oscillator (s), 

CAOs - creation and annihilation operators. 

CCRs - canonical commutation relations 

PM-operators - position and momentum operators 

QM - quantum mechanics; QFT - quantum field theory 

LA - Lie algebra, LAs - Lie algebras 

LS - Lie superalgebra, LSs - Lie superalgebras 

2. Lie (super)algebraic approach to quantum statistics 

We have already indicated that the most general approach to determine the admissible 
commutation relations between the PM-operators would be to find all common solutions 
of the compatibility equations (1.8). This task is however very difficult. For this reason in 
addition to the requirement the PM-operators to satisfy eqs. (1.8) we shall require these 
operators to generate a Lie superalgebra (LS) from the class A and more precisely the LS 
s/(3|iV). 

At this place one may ask why to introduce additional restrictions, postulating that 
the PM-operators generate a Lie superalgebra (LS)? And why a Lie superalgebra and not 
a Lie algebra or any other algebraic structure? One possible answer to this question would 
be to say that such an assumption was a good guess. And this would be not a wrong 
answer. There is however a deeper reason, which is based on two key observations (see 
also [33] for a more detailed discussions in the frame of both QM and QFT) . 

The first key observation belongs to Green. In 1953 he has shown that also the statis- 
tics of quantum field theory (QFT) can be generalized to what was later called parastatis- 
tics [34]. Green was also the first to realize that the infinitely many solutions found by 
Wigner in [5] are nothing but different inequivalent representations of one pair of para-Bose 
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(pB) creation and annihilation operators (CAOs) B"^ defined as 

B^ = ^{RTiP) ^ R=^{B+ + B-), P=^{B+-B-). (2.1) 

The generahzation for the Hamiltonian (1.1) goes as follows. Introduce in place of the 
PM-operators new unknown operators 

S± = (^) '^'kak T i(2m«a;n)- V^P^fe, (2.2) 

In terms of these operators the Hamiltonian and the compatibility conditions (1.8) read 
(see also [7]): 

^f=f EE{Si.B„-.}. (2.3) 

a=l i=l 

JV 3 

EE[i^£'^^.>'^i] = 2i?±, z = 1,2,3, a = 1,2,..., AT. (2.4) 

f3=l j=l 

Postulate that the operators B^- satisfy the triple relations 

[{Bl, B^}, S^fe] = SikSajie - OBI + SjkSpjis - rj)Bl). (2.5) 

It is straightforward to verify that the operators (2.5) satisfy the compatibility condition 
(2.4). Replacing the double indices with one index, ai — > /, etc one rewrites (2.5) in the 
form: 

[{Bj, S]}, = SiKie - OB'j + SjKie - r7)i?f ), ^,v,s = ±. (2.6) 

By definition the operators Bf are called para-Bose operators (pB-operators) . Hence 
these operators yield a new possible statistics for the oscillator under consideration. The 
para-Bose (pB) operators were introduced by Green in QFT for quantization of integer 
spin fields [34]. 

In the same paper Green [34, 35] generalized the Fermi statistics to para-Fermi (pF) 
statistics. The defining relations for any N pairs of para-Fermi CAOs read: 

[[Fh F]],Fk] = l{v- efSjKFl e)%KF], ^,v,e= ±. (2.7) 

Certainly the triple relations (2.6) and (2.7) are satisfied by Bose and Fermi operators, 
respectively. 
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The second key observation is that any n pairs of pF operators generate the orthogonal 
Lie algebra so{2n + 1) = S„ [36, 37]. In fact the linear span of F^, and [F^, F^], 
j,k = 1, ...,n, is already closed under further commutations as this is evident from (2.5). 

The LA Bn belongs to the class B of simple Lie algebras. There are four infinite 
classes of simple Lie algebras A,B,C,D (in order to avoid confusion we denoted them 
with italic letters). Therefore the Fermi and the para-Fermi statistics can be called B- 
statistics. In [38] it was shown that to each such class there correspond statistics (perhaps 
more than one): A—, S— , C— and D— statistics, which are appropriate for quantization of 
spinor fields (for A— statistics see [39]). 

Similarly, if one considers the pB operators as odd elements, then the linear span of 
all pB operators Sf, i = 1, ...,n, and all of their anticommutators close a Lie 

superalgebra [40], which is isomorphic to the basic Lie superalgebra osp{l/2n) = B{0\n) 
[41] in the classification of Kac [6] . 

The LS B{0/n) belongs to the class B of the basic Lie superalgebras [6]. Therefore 
the Bose and, more generally the pB statistics can be called B— statistics. Also in this case 
there exist four infinite classes of basic Lie superalgebras, denoted as A, B, C and X>, and 
again with each such class one can associate statistics. Since every Lie algebra is a Lie 
superalgebra (with no odd generators), the classes A, B, C, D are contained in A, B, C and 
T>, respectively. 

So far only the physical properties of ^ [1, 2, 7 - 10] and B-statistics [11, 12] were 
studied in more detail. Very recently however large classes of new solutions of compatibility 
conditions (1.8) for all basic LSs were determined [42] based on the results on generalized 
quantum statistic previously defined in [43]. 

Motivated by the above results, in the present paper we study the properties of a 
new class of Wigner oscillators, which can be called sl(3\N)-Wigner oscillators. This is to 
indicate that the position and the momentum operators of the particles are odd generators 
of the Lie superalgebra sl{3\N) (in appropriate (3 jA?")— modules) and generate it. But 
this will be the topic of the next section. 

3. S1(3|N) Oscillators. Representation independent results 

In this section we perform the main step towards explicit construction of si (3 1 A)- Wigner 
oscillators: we find common solutions of Eqs. (1.5), (1.6) and (1.8) with position and 
momentum operators which are odd elements in the LS sl{3\N) and generate it. 
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As a first step we replace the unknown operators Rak and Pak with new unknown 
operators, writing down the time dependence exphcitly: 



^/ \ / iV — 3|mQ,a; - , . . /I-^ — 31 . 

E{t)k,a+3 = \ TT R{t)ak - IS J ——r P{t)ak, 



N AN — ZrUoiOJ - , , /liV — 31 - , , 



R{t)ak — \ j I _ „| {E{t)k,a+3 + -E'(0a!+3,fe), 



777- coll/ 

P{t)ak = is J |^"_ g| (-E(t)A;,a+3 " E{t)a+3,k)- 



(b) Hamiltonian 

AT 3 

\N-3 



(3.1) 



where A; = 1, 2, 3, a = 1, 2, A?" and 

£ = +1 if iV>3 and £ = -1 if iV = 1, 2. (3.2) 
In terms of the new variables we have: 

(a) Position and momentum operators 



(3.3) 



N 3 

^ = E E nr^^{^Wfc,«+3, £^(t)a+3,fc}. (3.4) 



a=l k=l 

(c) Hamiltonian equations 

E{t)k,a+3 = i£'jjE{t)k,a+3: E{t)a+3,k = -i£<^-E'(^)a+3,fe, (3-5) 

(d) Heisenberg equations 

AT 3 

mk,a+3 = E ^[{£;(t),- ^+3, E{t)0+s,j}, E{t)k,a+3], (3.6a) 

N 3 

E{t)a+3,k = TTf- 

E E[-t^(^)^-.^+3' ^(^)«+3,fc], (3.66) 

(e) Compatibility conditions: 

^ JV 3 

E{t)k,a+3=j^J2J2^W)i^f^+^^^(^h+s,j}.E{t)k,a+3l A^3 (3.7a) 

^ N 3 

E{t)^+s,k - -^rr^^ E J2^{E{t)j,p+s, E{t)p+,j}, E{t)a+3,k], N ^ 3 (3.76) 



A-3 ^ 
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The time dependence of E{t)k,a+3 and E{t)a+3,k is evident from (3.5), despite of the fact 
that these operators are still unknown: 

E{t)k,a+3 = -E^fe,a+3(0)e'^'^^ £?Q,+3,fc(*) = -£^a+3,fe (0)e"~'^'^*. (3.8) 

From now on we set: 

-^'A;, 0+3(0) = Ea+3,k{Q) = -S'a+S.fe- (3.9) 

Then 



Rakit) = JrvrA, (Ek,a+3e'''^' + E^+s,ke-''^'), (3.10a) 

\/ \N — 3 mrvo; \ / 



\N — 3 mao; 



Pakit) = is J (Ek,a+3e'''^' - £;a+3,fee---*) , (3.106) 



whereas the Hamiltonian is time independent: 



N 3 ^ 



^ = X] 2Z liV - 31 Ea+3,k}- (3.11) 



a=l k=l 



Let us underline that the above equations (3.3)-(3.7) were obtained from (1.1), (1.5) 
and (1.6) just as a result of change of variables (3.1). E'a+a,^ and E'^^a+s are still unknown 
operators. As a first step now we find operators which satisfy Eqs. (3.7) and are odd 
generators of sl{2>\N). 

For convenience we consider sl{?>\N) as a subalgebra of the general linear LS gl{?>\N). 
The latter is a complex linear space with a basis Eab, — 1, 2, A^+3. The ^2-grading 
on gl{2>\N) is imposed from the requirement that 

EiAi Eai, ^ = 1,2, 3, A = 4, 5, A?" + 3 are odd generators, (3.12a) 
Eij, Eabi hj — ^,^,'^, A, B — 4,5, ...N + 3 are even generators. (3.126) 

The supercommutator on gl{3\N), turning it into a LS, is a linear extension of the relations 
{Eab, Ecd] = SbcEad - {-iy'''^^^^^'''''^^^°^SADEcB, (3.13) 

where 

{Eab, Ecd] = EabEcd - {-IY^^^^^''^''^^^^^''^EcdEab. (3.14) 
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In (3.13) and (3.14) A, B,C,D = 1, 2, N + 3. 
The LS sl{3\N) is a subalgebra of gli3\N): 

sl{3\N) = span(^gAEAA - 9bEbb, Ecd\C ^ D, A, B,C,D= 1, N + 3 



;), (3.15) 



where 



91 = 92 = 93 = 1, 94 = 95 = ■■■ = 9N+3 



= -1. 



(3.16) 



Then 



n' = spa.nl Eaa\A= l,...,iV + 3 



(3.17a) 



and 



H = span 



(^9aEaa-9bEbb\A,B=1,...,N + 3 



(3.176) 



are Cartan subalgebras of 9l{3\N) and sZ(3|A^'), respectively. These algebras are certainly 
commutative. 

It is evident from (3.12b) that the even subalgebra of gl{3\N) is the Lie algebra 
^/(3) © gl{N) with 

9I (3) = span{Eij\i,j = 1,2,3} and 9l{N) = span{EAB\A, B = 4,5, N + 3}. (3.18) 

Coming back to the Wigner problem we observe that the generators Ek^a+3, Ect+3,k 
with = 1,2,3, a = 1,...,N do satisfy the compatibility conditions (3.7). From now 
on we consider this particular solution, i.e., we identify the variables (3.7) with the odd 
generators (3.12a) of s/(3|A^). It is easy to verify that the linear span of E^^ce+s, E^+s^k and 
their anticommutators yield s/(3|A^). Hence, the position and the momentum operators are 
also odd generators (see (3.3)) and they generate s/(3|A'"). It is straightforward to verify 
that the time dependent operators (3.10) yield a simultaneous solution of the Hamilton's 
equations (1.5) and the Heisenberg equations (1.6). Hence Ri, . . . , Rjv and Pi, . . . , Pjv 
obey propositions (P5) and (P6). 

Already now we can draw certain conclusions which are consequences only of the 
postulates (P5) and (P6). These conclusions are representation independent, they have to 
hold in every state space. 

Using only the supercommutation relations (3.13) one derives from (3.11) 




JlUJ 



{NEii + NE22 + NE33 + 3E44 + ... + 3En+3,n+3) 



(3.19) 



N-3 
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and therefore the Hamiltonian is an element from the Cartan subalgebra Ti of sl{3\N). So 
are the Hamiltonians Ha of each individual particle, 

^a = ^+^R^ = ^^(£^ll + S22 + S33 + 3£;a+3,«+3), « = 1, iV, (3.20) 

and therefore the energy of each individual particle is preserved in time. Eq. (3.20) follows 
from 

- |^_3|^^^ (^M + ^a+3,a+3), (3-21) 
Pai = _ 3| + -E'a+3,a+3), (3.22) 

~ I AT- _ q| (-^11 + -^22 + -E'33 + 3-E'a+3,Q:+3)) (3.23) 

Pa = ^^(^11 + ^22 + £^33 + 3£;a+3,a+3). (3.24) 

In the above expressions i = 1,2,3 and a = 1,2, ...,7V. The conclusion is that indepen- 
dently on the representation and hence in every state space the set of all operators 

H, Ha, Rl, Pi, Rli, Pli, z = 1,2,3, a = 1,2,..., AT, (3.25) 

constitute a commutative set of operators, whereas neither coordinates nor momenta com- 
mute, see (1.2). In this respect the WQO belongs to the class of models of non-commutative 
quantum oscillators [44 - 48] and, more generally, to theories with non-commutative geom- 
etry [49], [50]. Moreover, as in [51], the coordinates of the particles are observables with a 
quantized spectrum just like energy, angular momentum, etc. 

This is an appropriate place to mention that any physical observable, which is a 
function of only even generators is time independent, which is another representation 
independent property. The latter stems from the observation that the even generators 
commute with the Hamiltonian, 

[H,Eij\=Q for any i,j = 1,2,3 or i, j = 4, 5, TV. (3.26) 

An important physical observable of this kind is the angular momentum. For the 
components of the orbital momentum M^i of each individual particle a = 1,2, N, we 
postulate the same expression as in the canonical case: 

1 ^ 

~ 2 ^ ^ikl{Rak, Pal}- (3.27) 
k,l=l 
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Being an anticommutator of odd observables Rai and Pqi) each M^i is an even element and 
therefore it is an integral of motion, it does not depend on the time t. It is far from evident 
however that the above definition is physically acceptable. One necessary requirement is 
that for each a Mai, Ma2, transform as vector operators under space rotations. The 
position and the momentum operators of each particle have to transform also as vectors 
with respect to space rotations. This is clear. But so far we do not have generators of 
the physical rotation group. In order to determine them we proceed also as in canonical 
quantum mechanics. 

To begin with we express the projections of the angular momentum of each particle 
via Weyl generators: 

eh ^ ^3 
Mai = XI SikiEki = -i ^_g X SikiEki. (3.28) 

' ' k,l=l k,l=l 



Setting 



we obtain: 



Mai = Ji^Sai (3.29) 



Sai = — i X £iklEkl- (3.30) 



Explicitly 



'S'al = i(-E'32 — -E'23), 'S'q,2 = i(-E'l3 — -E'31), 'S'a3 = ^(-E'21 — -E'12) (3.31) 

with commutation relations 

3 

[SgjjSgk] — ^'y^^SjklSgl, (3.32) 

/ = 1 

which are the known commutation relations between the components of the angular mo- 
mentum also in conventional QM. There is however one essential difference. In the canon- 
ical case the operators Sai, 3^2, S^s, measure the components of the angular momentum 
of any particle in units h. Here they are measured in units h/\N — 3\. 

The most striking difference between the canonical oscillator (or the sl{l\3N) oscillator 
[9, 10]) and the s/(3|A^) oscillator comes from the observation that the angular momentum 
operators Sgi, Sa2i do not depend on a. They are the same for all N particles (see 
the RHS of (3.30), (3.31)). In the next sections we will discuss this feature in more detail. 
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For the components of the angular momentum of the oscillator we have 



^ "'^ N-3 

a=l k,l=l 



where 



or 



Sj = -i ^ EjkiEki, (3.34) 

k,l=l 

Si — i{Es2 — E23) , S'2 = i(-E'i3 — -Esi), = i(£'2i — -E12). (3.35) 

In view of (3.33) 5'^ measures the components of the total angular momentum of the 
oscillator in units hN/{N — 3): 

3 

[Sj,Sk]^i'YejkiSi. (3.36) 



1^1 

The operators 

S+ = Si + iS2 = i-E'32 — i-E'23 — Eis + i?3i, (3.37) 

S- = Si — 15*2 = i£'32 — i-E'23 + £"13 — £"31, (3.38) 

-53 = i(£2i-£i2), (3.39) 

satisfy the known commutation relations for the generators of the algebra so(3). 

[4, S+] = S+, [Ss, 5_] = [S+, 5_] = 2^3. (3.40) 

At this place we postulate that Si, S2, Ss are the generators of space rotations. One verifies 
that 

333 

[Sj, Mak] = i'Y^^jklMal, [Sj, Rak] — i'Y^^jklRal, [Sj,Pak]—i'Y^^jklPal, (3.41) 
1=1 1=1 1=1 

i.e., the components of M^, Rq,, Pq, of each particle transform as vector operators with 
respect to space rotations. This holds in any representation. Moreover the total Hamil- 
tonian and the Hamiltonians of each individual particle are scalars with respect to space 
rotations: 

[Sj,H] = [Sj,Ha] = 0, j = 1,2,3, a =1,2,..., AT. (3.42) 
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From the results obtained so far we can draw some further conclusions. First of all 

[H, Pi] = [H, RU = 0. (3.43) 

The eigenvalues of should be interpreted as squares of the admissible values for the ith 
coordinate of particle No a. The circumstance that -R^j commutes with the Hamilonian 
then means that the square of the ith coordinate of the ath particle is an integral of 
motion. Since, moreover, all operators and commute with each other 

[Pai^P^j] = [PaiiR^fij] = [P-aiiR^fij] = 0) (3.44) 

they can be measured simultaneously. Observe that the above statement is representation 
independent, it has to hold within every admissible state space. 

4. S1(3|N)-Wigner quantum oscillators. State spaces 

So far we have introduced time dependent operators Ri, . . . , Rjv and Pi, . . . , Pjv which 
obey postulates (P5) and (P6). In the present section we determine state spaces, which are 
simultaneously representation spaces of sl{3\N), so that the oscillator becomes a WQO. 
In principle one could search among all representations of sl{3\N) and select those of 
them for which the Wigner problem has solutions. Since however explicit expressions 
for all representations are not available, we restrict our considerations to the class of 
ladder representations of s/(3|A^) [52] (leaving for future another class of known irreps, the 
essentially typical representations of s/(3|A^) [53]). 

Below, following [52], we recall shortly the main properties of the ladder representa- 
tions directly for sl{3\N). Let 

1. {cf = bf, = 6^, = 6^} be Bose operators considered as odd elements: 
[K,b+] = Sij, [b+,b+] = 0, [br,bj] = 0, i,j = 1,2,3; 

2. {cf = f^, cf = f^, . . . , = /m+s) t»e Fermi operators considered as 
even elements: 

{frJj'} = Sij^ = 0, {/-,/-} = 0, z,j = 4,...,Ar + 3 (4.1b); 
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(4.1a) 



3. 



The Bose operators anticommute with Fermi operators. 



(4.1c) 



For a justification of this unusual grading see [54]. Here we only remark that with 
this grading the Fock space introduced below gives an infinite-dimensional irreducible 
representation of the orthosymplectic LS B{m\n) = osp{2m + l|2n) with generators the 
Bose and the Fermi operators and their supercommutators according to the grading. 

Denote by E{?>\N) the Bose-Fermi algebra, namely the free superalgebra, generated 
by the CAOs c^, . . . , c^_|_3 with the relations (4.1). 

We are going to work in the Fock module W{2>\N) of E{2>\N) with an orthonormed 
basis 

\n„n„...,n^+s)= ^'^ ' 7 10), c^|0) = 0. (4.2) 

Vni!n2!n3! 

where 

ni,n2,n3eZ+, 77,4, 77,5, tt-at+s e {0, 1}. (4.3) 

For definiteness ni,n2,n3 are said to be bosonic coordinates of the state \p',n), whereas 
77-4, 77-JV+3 are referred to as fermionic coordinates of \p; n). We shall see that each basis 
vector \p; n) determines up to a sign the possible coordinates of the particles. We call the 
basis (4.2) a Fock basis. 

Clearly the representation of E{S\N) in W{S\N) is infinite dimensional. The trans- 
formations of the basis (4.2) under the action of the CAOs cf read ([52], Eq. (48)): 



6+|..,ni, ..) = Vni + l\..,ni + l,..), 1,2,3; (4.4a) 

b~\..,ni,..) = ^/f^i\..,ni-l,..), i= 1,2,3; (4.46) 

/+|..,n„ ..) = (-l)"i+-+"*-ix/r^|.., m + 1, ..), i = 4, 5, .., N + 3; (4.4c) 

f-\..,ni, ..) = (-l)-^+--+"*-i V^|..,ni - 1, ..), z = 4, 5, ..,N+3; (4.4d) 

19 



It follows ([52], (49), (50)): 



c+cj\..,ni, .., nj, ..) = (^i)"^+-+'^-i (^,)-i+-+"^-i -^(1 + gini)nj\.., m + 1, .., n,- - 1, ..), 
for z < j = 1,2, ...,Ar+3; (4.5a) 



c+c-\..,nj,..,ni,..) = {g,r^+--+^^-^-\gjr^+--+^^-^^^^ 

for i>j = 1,2, ...,7V + 3; (4.56) 

c+c7|..,nj, ..) = nj|..,ni, ..), i = 1, 2, A/" + 3. (4.5c) 
Explicitly Eqs. (4.5) read: 

b+bj\..,ni, ..,nj, ..) = y(n~Tl)n~|..ni + 1, ..,nj - 1,..), i <j = 1,2,3; (4.6a) 

bfbj\..,nj, ..,ni, ..) = y^(ni + l)nj|..nj- - 1, ..,ni + 1, ..), i> j = 1,2,3; (4.66) 

/+/r |.., m, .., n,-, ..) = (-l)n.+■.+n,-^ ^(1 _ n,)n,- 1 ..n^ + 1, .., n,- - 1, ..), 

j >z = 4,5, ...,Ar+3; (4.6c) 

/+/- n,-, .., n,, ..) = (-l)".+-+-^-i-y (1 - ni)nj\..nj - 1, .., n, + 1, ..), 

i > j = 4,5, ...,iV + 3; (4.6d) 

6+/r|.., m, .., n,-, ..) = (_i)-i+-+-.-y(^. + i)^.|..„. + 1, ^. _ 1, ..), 

z = 1,2,3, j = 4,5, ...,iV + 3; (4.6e) 

f+b-\.., nj, .., n,, ..) = (-l)-i+-+-*-y (1 - nOn,-|..n,- - 1, .., + 1, ..), 

z = 4,5...,iV + 3, J = 1,2, 3. (4.6/) 

It is straightforward to verify that the following proposition holds [52] . 

Proposition 4.1: The linear map defined on the generators as 

Eij^c+cJ, i,j = l,2,...,iV + 3, (4.7) 

gives a representation of the LS gl{3\N) in W{3\N). 

From now on we write Eij also for cfcj, i.e., we use for simplicity one and the same symbol 
for the abstract generators of gl{3\N) and for their images as operators in W{3\N). 
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From Eqs. (4.5) one concludes that the infinite-dimensional Fock space 1^(3|A'") re- 
solves into an infinite direct sum 

oo 

Wi3\N) = ^ViN,p), (4.8) 

p=0 

of finite-dimensional subspaces 

V{N,p) = span{\ni, ...,nN+3) \ rii + ... + njv+3 =p}, (4.9) 
labelled with p e N (all positive integers). We set 

|ni, n2, ns, n^, un+s) = \p; ni,n2, ns, n^, un+s), (4.10) 

if we wish to underline that |ni, ...,njv+3) £ V{N,p). 

Each subspace V{N,p) is invariant (and in fact irreducible) with respect to the opera- 
tors (4.5) and hence with respect to any physical observable. Therefore each such subspace 
is a candidate for a state space of the system. 

One verifies that the subspace 

V{N,p,nb,n^,...,nN+3) cV{N,p), (4.11) 

which is a linear span of all vectors \ni, ...jUn+s) with ni + 77-2 + ns = rib being fixed and 
n4, ...,njv+3 being also fixed is an irreducible gl{3) module. Similarly, the subspace 

V{N,p,nf,ni,n2,n3) C V{N,p), (4.12) 

which is a linear span of all vectors |ni, un+s) with both ni, 77-2, and n/ = 77,4 + ... + 
njv_(_3 fixed is an irreducible gl{N) module. Finally, the subspace 

V{N,p,nb,nf)GV{N,p), (4.13) 

which is a linear span of all vectors \p',ni, ...jUn+s) with both 77.1+77-2 + ^3 = rib and 
71/ = 7i4 + ... + Tijv+s is an irreducible gl{3) ® gl{N) module. The labels ni, rif and p in 
(4.13) are not independent since rib + rif = P- Nevertheless we prefer to keep the more 
symmetrical notation in (4.13). 

Proposition 4.2. The N— particle 3D oscillator with PM-operators (3.10) and a state 
space V{N,p) can be turned into a Wigner quantum oscillator for any p = 1,2, .... 
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Proof. 



(1) . Every finite-dimensional linear space with a scalar product and in particular V{N,p) is 
a Hilbert space. Postulating that to every state of the system there corresponds a normed 
to 1 vector from V{N,p) one fulfills the first requirements (PI) of the definition of a WQS. 

(2) . Eqs (4.4) yield that the Hermitian conjugate of a creation operator is an annihilation 
operator: 

{b+r^b-, ^ = 1,2,3, {fXr-fl, A = 4,5,...,iV+3. (4.14) 

Therefore 

E*^ = ictc-r = ic+c-) = Ej,, i,j = l,...,N + 3. (4.15) 

From (3.10) and (4.15) one concludes that the PM-operators are Hermitian operators in 
V{N,p). As a consequence also the Hamiltonian (3.11), the Hamiltonians of each individual 
particle, the projections of the angular momentum M^j and Mj are Hermitian operators. 
Hence the requirement (P2) holds too. 

(3) The validity of (P5) and (P6) was already established in the previous section. 

(4) Finally we postulate that any observable L can take only values which are eigenvalues 
of L (P3) and that the expectation value of L in a state ip is evaluated according to (P4). 

This completes the proof. 

In this way to every p — 1,2, ... there corresponds an sl{3\N) Wigner quantum os- 
cillator with a state space V{N,p). The PM-operators corresponding to different p are 
inequivalent because they correspond to different irreducible representations of sl{3\N). 

5. Physical properties - energy spectrum 

In this section we begin to discuss the physical properties of the sZ(3|A^) WQOs. We 
compute the energy spectrum of the system and of the individual oscillating particles. We 
shall see that the energy spectrum is equidistant, but finite. Another surprise is that in 
certain state spaces each individual particle can have a zero energy. 

Let V{N,p) be ap-state space, see (4.9). All vectors |p; ni, n2, na; 71,4, ...,71^+3) (with 
ni + ... + un+s = p) constitute a basis of eigenvectors of H. For the eigenvalues of the 
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Hamiltonian (1.1) in this state one obtains from (3.19) and (4.5): 

3 Ar+3 



Therefore the energy £'(p; ni, nAr+s) of this state is 

3 7V+3 

|iV-3 



E{p- m, nAT+g) = (A^ J] + 3 J] n^) . (5.2) 



Clearly all vectors |p; n) = \p;ni, un+s) with one and the same 

nb = ni + n2 + ns and Uf = + ... + un+s, (5.3) 
have one and the same energy: 

E{N,p,nb,nf) = j^^^^{Nnb + 3nf). (5.4) 

The energy E{N,p,nb,nf) depends actually on three independent variables, for instance 
N,p,nf since rib = p — Uf. 

Replace in the RHS of (5.4) Ub with p — nf. 

E{N, p, Ub, Uf) = {Np - (AT - 3)nf) = uh - ^^n/) . (5.5) 
Taking into account that 

n/ = 0, 1, 2, min(A'",p) <^=^ Ub = p,p — 1, max(0,p — N), (5-6) 
one concludes: 

Corollary 5.1 The spectrum of the Hamiltonian H in V{N,p) consists of 

d{N, p) = mm{N, p) + 1 (5.7) 

equally spaced energy levels, with spacing u>h. The energy E{N,p,nb,nf) corresponding to 
a particular value of rif is (5.5). 

Here comes the first big difference with a system of A?"— particle free canonical 3D 
oscillators. In the latter case the energy is also equally spaced with the same spacing %uj. 
Now however there is an infinite number of energy levels: 

Eq = nuj{^N + (^, 5 = 0,1,2,.... (5.8) 
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In order to determine the multiplicity of the energy E{N,p, n^jUf) we have to compute 
the dimension of the eigenspace V{N,p,nb,nf) of H corresponding to this eigenvalue. 
By definition, see (4.13), 

V{N,p,nb,nf) = span{\p;ni, ...,nN+3) \ ni + n2 + ns = Ub, Ub + Uf = p}. (5.9) 

All states from V{N,p,nbjnf) are stationary states, they have one and the same energy. 
Each subspace V{N,p,nb,nf) is a gl{3)(Bgl{N) module. There are d{N,p) — min{N,p) + l 
such modules, 

min{N,p) 

V{N,p)= J2 ®V{N,p,nb,nf), nb=p-nf. (5.10) 

nf=0 

With respect to the even subalgebra each V{N,p,nb.,nf) behaves as a tensor product 

V{N,p,nb,nf) ^Vi{N,p,nb)(^V2{N,p,nf), Ub + Uf^p (5.11) 

of a (3)— module Vi{N,p,nb) and a ^Z( A?")— module V2{N,p,nf). 

The linear space Vi {N, p, Ub) is a Fock space of three Bose operators with a basis 

ic^^^^ ic^^^^ ic^^'^^ 

|ni,n2,n3) = ^^' y ^ , — |Q), c^|0) = 0, m + n2 + n^ = nb = p-nf. (5.12) 

We say that Vi(A", p, n^) is the bosonic component of V{N,p, Ub, Uf) or a bosonic subspace. 
It is a simple exercise to verify that Vi{N,p,nb) is an irreducible gl {3)-module: there is 
only one eigenvector of the Cartan subalgebra of gl{3) annihilated by the positive root 
vectors. This vector is the highest weight vector |n6,0,0) with a highest weight {ub, 0,0). 
The dimension of Vi {N, p, Ub) is 

dimVi(iV,p, Hb) = {nb + l){nb + 2)/2. (5.13) 

Clearly the algebra so(3) of the rotation group SO{2>) is a subalgebra of the algebra 
gl{2>) (see the expressions for the so(3) generators(3.37) - (3.39)). Therefore the space 
rotations transform only the bosonic part yi(A", p, rib) of V{N,p, nb,nf). 

The fermionic subspace V2{N,p,nf) (the fermionic part of V{N,p, Ub, Uf)) has a basis 

(/4)"^(/5)"^-(/iv+3)""+^|0), n4, ns, . . . , e {0, 1}, 71^ + 71^ + ... + riN+s = Uf. 

It is also irreducible gl{N)— module with a highest weight vector |li, I2, ...In/, 0..., 0) and 
dimension 

AT! 

dimV'.(iV,p,„,) = -j^^^-_^. (5.14) 
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Therefore 

dimV(iV,p.„.,„,)=(!^x-^^^. (5.15) 

Observe the pecuharity of only one ?>D Wigner osciUator (A^ = 1): for any p the 
oscillator has only two energy levels, namely E = ujTi^ and E = ijjh{'^ + 1). In the case 
p=l the result reduces to the first two energy levels of a canonical 3Z) oscillator. 

Let us summarize. 

Corollary 5.2 The linearly independent states in V{N,p) corresponding to the energy 
E{N,p,ni),nf) are all vectors |p; ni, niv+a) with 

3 JV+3 

^ni=nb, ^nA = nf (5.16) 

i=l A=4: 

(and rib + nf =p). Their number, and hence the multiplicity of E{N,p,ni„nf) is given by 
Eq. (5.15). 

For the ground energy E[N.,p)min of the system one derives: 

If iV > 3, iV > E{N,p)^in = ]^3p, m = ^^"l^") 
If iV > 3, AT < p, E{N,p)rmn = J^N{p - iV + 3), m= (5-176) 
If AT < 3, E{N,p)min = ^^Np, m=^(p + 2)(p+l), (5.17c) 

where m denotes the multiplicity of E{N,p)min, the number of the linearly independent 
ground states. From the above results one concludes: 

Corollary 5.3 The ground energy of the system is always positive. In the cases N > 3 
the ground state is nondegenerate only if N = p. 

The circumstance that the energy of the ground states is never zero is not surprising. 
The same holds also in the canonical case. In conventional quantum mechanics also the 
energy of the ground state of each individual oscillator is positive. Is this the case for the 
WQOs? No, it is not. We proceed to show this. 

Since the Hamiltonians Hi,...,Hn of the oscillating particles are elements from the 
Cartan subalgebra of sZ(3|A'"), they commute with each other and with the Hamiltonian 
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of the system (1.1). Hence the energy of each individual particle is preserved in time. 
The energy of the a— th particle when the system is in the state |p; ni, tt-at+s) is the 
eigenvalues of Ha on this state. Since 

-^ab;"-l, AT+s) = i^^^-^ ("-1+^2 + 77-3 + 3na+3)b;ni,...,njv+3), (5.18) 

the energy of the a — th particle in the state \p;ni, njv+s) is 

Ea{p;ni,...,nN+3) = p^32|(m + ?72 + ?73 + 3na+3). (5.19) 

Observe that in a given state \p; n) all particles can have at most two different energies: 
all particles No CKjjCKj,... with fermionic numbers tt-q^ — ^oij — ••• 

= have one and the 

same energy £'0 = whereas the energy of the rest (those with fermionic numbers 

one) is El = ^|^(n6 + 3). 

As for the energy spectrum of each particle (after some combinatorics) one obtains: 

Corollary 5.4. The spectrum of the Hamiltonian in V{N,p), measured in units 
?kij/\N — 3\, or, which is the same, the spectrum of the operator Eii + E22+E3s + 3Eoi+3,a+3 
(see (3.20)), reads: 

{p-min{N,p),p-min{N,p) + l,...,p,p+l,p + 2}, iV > 3, p > 2; (5.20a) 
{0,1,3}, for AT > 2,iV7^ 3, p= 1; (5.206) 

{p- l,p,p+ l,p + 2}, forA = 2, p>2; (5.20c) 
{p,p + 2}, foriV = l,p>l; (5.20cZ). 



In all cases with p > 1 and N > 1 there is a finite number of equally spaced energy 
levels, with spacing one (in units j^^^)- Peculiarities appear however in the state spaces 
with p = 1 or N = 1. In the case (5.20d) the spacing is twice bigger compared to p > 1 
cases, whereas in the cases (5.20b) the equally spacing rule is violated. 



6. Physical properties - oscillator configurations 

Here we analyze the 3D space structure of the Fock states \p; n) . Our results are 
based essentially on postulates (P1)-(P6). We begin however with a few observations of a 
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more general character. In particular we indicate that the superposition principle holds for 
any WQO. We point out also that the dispersion (or the square root of it - the standard 
deviation) is a convenient tool to search for particular states where the observables can 
have simultaneously particular values. 

We are not going to discuss the details of the quantum measuring process and of how 
to prepare a system in a particular quantum state. Until now these are topics of increas- 
ingly hot discussions (see, for instance [55] and the references therein). We find however 
appropriate to recall the experimental definitions of certain entities and in particular of 
such in principle well known concepts as a mean value and a dispersion based on Gibbs 
ensemble (gedanken) experiments and to relate them with the theoretical predictions. 

The definition to follow is not the most general one. It is adjusted directly for our 
considerations. Let be a particular state of the quantum system under consideration. 
Then the Gibbs ensemble (GE) corresponding to this particular state consists of a large 
number A^o of identical quantum systems 

all of them prepared to be in the state ^. By ^^'^^ we denote the kth individual such 
quantum system. 

Consider the observable L. Perform with all quantum systems from the GE simulta- 
neously one and the same experiment, namely measure the value of the observable L in 
the state ^. Let ni quantum systems among all A^o registered a value Zi, n2 - a value I2, 
and so on, '^f.Uk = Nq. Then the average (=the mean = the expectation) value of L in 
the state ^ is 

{L)l^P = J2^l, = Y^PJ,, (6.2) 

where Pk = fik/No is the probability to measure the result Ik- Note that by construction 
^1 7^ ^2 7^ ••• It is assumed that the object under consideration admits a statistical descrip- 
tion. This means that the values of the probabilities Pk are practically independent on the 
number Nq of the identical systems in the ensemble if this number is sufficiently large. 

The first conclusion, which is due to postulate (P3), is that all numbers Zi, /2, •• have 
to be eigenvalues of the operator L. Next, the "experimental" expectation value of L (6.2) 
should be consistent also with postulate (P4). Let \E'g be an orthonormed system 

of eigenstates of L: L'^k — h'^k, k — 1, 2... so that 

* = + ^2*2 + + (6.3a) 
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Then postulate (P4) yields: 

^lyheor ^ ^^ l^^ ^ J2 \ak\%, (6.36) 

k=l 

This is actually the superposition principle for WQOs. 

Superposition principle: Let L be an observable and ^ be a normed to one state, which 
is a linear combination \1/ = ai^i+a2^2 + ■■■ + otg'^ of an orthonormed set of eigenstates 
*i,*2,---,*3 o/^- ^^'fe = /fe^'fe. Then 

where each coefficient \ak\^ = Pk gives the probability of measuring the eigenvalue Ik of L 
corresponding to the eigenstate ^k- 

A comparison of (6.2) with (6.4) gives that 

{L)%^^ = {1)^'°' = (L)^, (6.5) 

if Pk is a sum of all \a1\ for which 1^ = lk- In particular if the spectrum of L is nondegen- 
erate (h h •••)) ^^en |q;|| = Pk- 

If ^ is an eigenstate of L, = l'^, then the expectation value of L in the state 
yields the eigenvalue I, (L)* = L The inverse is in general not true: from (L)* = I it 
does not follow that L'^ = l'^. In other words from (L)^ = I one cannot conclude that 
the measured value of L in each single experiment is I. The entity which insures that all 
individual experiments measure one and the same eigenvalue for L in the state is the 
dispersion. The dispersion Disp(L)^ of L in the state ^ is by definition 

Disp(L)* = J^Pkih - {Lhf (6.6) 

k 

Then one verifies, taking into account the relation \ak\'^ — Pk, that 
g 

Disp(L)* = J2 Pkih - {L)^f = (L^)^ - (L)| = (*, L^*) - (*, L^f. (6.7) 

k=l 

As a consequence of (6.2) and (6.3) one has: 
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Proposition 6.1:The next three statements are equivalent: 
(a) ^ is an eigenstate of the observable L: = l^, 

i.e. the observable L has a definite (= a particular) value I in the state ^ . 
(h) The dispersion of the observable L in the quantum state ^ is zero, 
(c) All Nq individual experiments from the GE register one and the same eigenvalue for L. 

The generalization of proposition 6.1 to the case of any number of observebles requires 
some care. 

Proposition 6.2. The next three statements are equivalent: 

(a) ^ is simultaneously an eigenstate of each observables Li, L2, . . . , Lm-' = 1^^. 

(h) The dispersion of all Li, . . . , Lm in the quantum state ^ is zero. 

(c) In the state all observables Li, . . . , Lm can be measured simultaneously in each 

individual experiment from the GE. In each such experiment they register one and 
the same eigenvalue 1^ for each L^. 

The equivalence of (a) and (b) is evident from the very definitions (6.3a) and (6.6). 
To prove the part (c) it suffices to show that the measurement of any observable Lk does 
not disturb the simultaneous measurements of the rest of them, so that one can apply 
proposition 6.1 to any L^ separately. To this end we follow the considerations of Dirac 
[3] in QM. Without specifying the domain of definition of the operators, Dirac assumes 
that the simultaneous measurements of Li, L2, ■ ■ ■ , Lm do not disturb each other if these 
operators mutually commute. 

Let D G W he the linear span of all common eigenvectors (fi, ...,(/?^, of Li, ...,Lm- 
The subspace D is invariant with respect to the action of Li, ...,Lm- Moreover the oper- 
ators Li, ...,Lm, commute on D. Therefore, following the arguments of Dirac, we accept 
that the measurement of any observable Lf., whenever the oscillator is in a state from 
D, does not disturb the simultaneous measurements of the other observables. In partic- 
ular this holds for any eigenvector (pk- Therefore one can apply proposition 6.1 for any 
observable Lk separately, thus proving proposition 6.2. 

Coming back to the space structure of the oscillator, we recall that we work in a 
rectangular coordinate system. By 

e= (61,62,63) (6.8) 
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we denote the three orthonormed frame vectors. Unless otherwise stated by coordinates 
of any 3D vector we have always in mind the coordinates with respect to this frame. 

Our main goal in this section is to show that typically a state \p; n) can be interpreted 
as having a 3D-configuration with the property that one or more of the particles are allowed 
to occupy only a finite number of points. 

For convenience we shall work not with the operators Rak{t) and Pak{t) themselves 
but with their dimensionless version 

Takit) = \n^^4^^-k{t) = Ek,a+3e''^' + ^a+3,fee---*, (6.9a) 



Pak{t) = y = ie(^Ek,a+3e'''^' - E^+s^ke-''"^') . (6.96) 

The peculiarity of the WQOs stems from the observation that the geometry of these 
oscillators is noncommutative: 

[rai,r/3j] 7^ 0, ipai.Pfdj] 7^ 0, (6.10) 

whereas in QM they do commute. On the other hand however all 6A'" operators f^j, p^j, 
i,j = 1, 2, 3, a,P = 1, 2, N, do commute with each other: 

[rli,r^pj] = [fli,Ppj] = \pli,ph] = 0, (6.11) 

whereas in QM they do not commute and more precisely ^ 0. For this reason, as 

already mentioned in the introduction, we say that the geometry of the Wigner oscillator 
is noncommutative but square commutative. 

Proposition 6.3. The operators 

rl, = Eu + E^+s,a+3, ^ = 1,2,3, a = 1, 2, iV, (6.12) 

constitute a complete set of commuting operators in V{Njp). In particular the eigenvalues 
r^j of f'^^ determine uniquely the basis vectors |p; ni, n2, ns; ^4, n^v+s) 

Proof. A given subspace V{N,p) consists of the linear span of all vectors 

|p;ni, ...,njv+3) with m + n2 + ns + n4 + ... + nN+3 = P- (6.13) 
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Any vector (6.13) is an eigenvector of f^^: 

^ai |P; "-1 ' • • • ' "-iv+s) = {rii + ria+s) b; ni , . . . , njv+3) • (6.14) 
The claim is that all eigenvalues 

rli = rii + na+3, i = 1,2,3, a = l,2, ...,iV (6.15) 

of f^^ determine uniquely the vector \p;ni, ...jUn+s)- 
Set 

N+3 

E ''Ai = ai. (6.16) 

Then 

N+3 

r^j = iVrii + n4 + ns + ... + riiv+s = Oi, i = 1, 2, 3. 

In order to eliminate 77.4, nAr+3 add in both sides of (6.16) ni + n2 + n^: 

Nui + ni + n2 + ... + nN+3 = flj + ni + n2 + na, i = 1, 2, 3. 

Taking into account (6.13) we obtain three equations for three unknown tt-i, 77-2, entities: 

(TV - l)ni - n2 - n3 = ai - p, 
-ni + (N - 1)77-2 - nz = a2 - p, 
-ni -n2 + {N - 1)773 = 02 - 



Their solution reads 



{N -2)ai + a2+a3-pN 
ni = 



N{N - 3) 
{ai + {N -2)a2 + a3-pN 
= NiN^) ' ^^-'^^ 

_ ai + (12 + {N - 2)a:i - pN 
~ N{N - 3) 

The values for the other entities Ua+s, a — 1,2, N follow from (6.15): 

na+3 = ria-ni, i = 1,2,3, a = l,2,...,N. (6.18) 
This completes the proof. 
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We proceed next to study the space configuration of the system whenever it is in a 
basis state \p; n) . Any such state is simultaneously an eigenstate of all operators f^j = 

Pai = ^ii + Ea+3,a+3' 

rli\p;-,ni,...,na+3,-) = r^Jp; .., n^, na+3, ••) (6.19) 

and of the Hamiltonian. Therefore according to proposition 6.2 all these 6A^ operators can 
be measured simultaneously in each individual quantum system from the ensemble. The 
eigenvalue r^^ of f^^ on \p; n) yields the square of the ith coordinate of particle a: 

rli ^ rii + na+3, a = 1, 2, .., A^, i = 2; (or 1, 2, 3). (6.20) 

From the last resut (6.20) one concludes. 

Corollary 6.1. A state \p;n) corresponds to a configuration when simultaneously for 
k = 1, 2, 3 the k — th coordinate of ath particle is either y/rik + ria+a or — -^nfc + ria+s- 

What (6.20) does not say however is what is the probability the fcth coordinate to be 
^/nk + rioi+s or —^rik + ?^a+3. The next proposition answers this question too. 

Proposition 6.4. // the system is in the state \p;n), then with equal probability 1/2 
the first coordinate of particle a is (measured to be) either y/ui + n^+s or —\Jn\ + na+s, 
the second coordinate of the same particle is either -\/n2~-Fri^+3 or — y'n2 + tIq,+3 and 
the third coordinate is either y^ns-Fn^H-s ^'^ — -\/n3 + n^+s. Also with probability 1/2 the 
kth component of the momentum of particle a take values ^/uk + Ua+s or —^/uk + Ua+s, 
k = 1,2,3.. 

Proof. The proof is based on the superposition principle and the explicit expressions for 
the eigenvectors of the coordinate operator fa,k, k = 1,2, 3. The latter read: 

a. Eigenvalue of fctk '■ 0. 

Eigenstates : 'i;°fc(.., Ofe, .., 0^+3,..) = |p; .., 0^, 0^+3, ..), (6.21a) 

b. Eigenvalues of rak '■ ±V^fc (^fc 7^ 0) : 

Eigenstates : v^,^{.., n^, .., Oa,+3,..) = (\p; ■■, Uk, .., 0^+3, ..) 

T (-l)"i+-+-«+2e-'='^*b; .., Uk - 1, .., la+3, ..), rik > 0, (6.216) 
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where in place of the unwritten indices one inserts all admissible indices which are the 
same in the RHS and the LHS of one and the same equality. The vectors (6.21) constitute 
an orthonormed basis of eigenvectors of rak{t) in ViN,p) for any a and any k. 
The inverse to (6.21b) relations take the form: 

|p;..,n,,,..,n«+3,..) = -L(-l)K+-+-«+^+i)-«+3eisn„+3a.t 

(^afc("' + na+3, 0^+3, ••) + {-1^"+' V+^{.., Uk + n«+3, 0^+3, ..)), (6.22) 

Note that (for any k = 1,2,3) the absolute values of the coefficients in front of and 
v~j^ in (6.22) are equal and their square is 1/2. Then the superposition principle asserts 
that if ffc, = 1, 2, 3 is measured, then with equal probability 1/2 the ctth particle will be 
found to have a kth coordinate v^nfe~-Fn^+3 or —^/nk'+n^^, respectively. 

The proof for the momentum is similar (see appendix E). This completes the proof. 

Let us underline. For a given state \p;n) from V{N,p) the interpretation of 

rak = ±^/nk + na+3, k = 1,2,3, (6.23) 

as coordinates of the particle a make sense only because |rQ,i|, |rQ,2| and Ir^sl are measured 
simultaneously in every individual experiment from the Gibbs ensemble. The circumstance 
that the coordinate operators fai, and f^s do not commute and therefore cannot 
have particular values is not used for the conclusion. The state \p; n) is anyhow not 
(and cannot be) an eigenstate of the coordinate operators. The conclusions are based 
on the fact that \p; n) is an eigenstate of the squares of the coordinate operators. The 
probability distribution for the coordinates, based on the superposition principle also does 
not contradict to the conclusions made so far. 
We denote as 

T{\p;n),a) = {±^^ni + na+3 ei ± y^n2 + na+3 e2 ± y^n^ + na+a es} (6.24) 

the positions where the ath particle can be measured to be (we say also "where the cuth 
particle can be accommodated" ) . 

With equal probability 1/2 the ath particle will be measured to have a A;— th coordinate 
^/nk'+n^^ or —^Jn]^ + nQ.+3, k= 1, 2, 3. Therefore the particle cannot be localized in only 
one of the points (6.23). 
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Following [9] and [10] we call the positions (6.24) nests. One should remember that 
the nests are not elements from the state space. They are just places where the particle 
can be measured to be. As we shall see, for certain states \p; n) the nests (6.24) do not 
exhaust all possible nests. In such a case the probabilities in proposition 6.4 are conditional 
probabilities. 

In the remainder of this section we will concentrate mainly on the determination of 
all nests corresponding to a given state \p; n) . On the way we derive an analogue of the 
uncertainty relations. We begin with an example. 

Example 6.1. Let N=l. Consider a state (f = \p = 12; 3, 0, 8; 1) and let the experiment 
measures simultaneously the absolute values of the coordinates x, y, z of the nests of the 
particle. Then each experiment from the Gibbs ensemble gives |ri| = 2, \r2 \ = 1, |r3| = 3, 
which means that ri = ±2, r2 = ±1, = ±3. Hence there are 8 nests where the first 
particle can be accommodated: 

r(|p; 3, 0, 8; 1, O^+g), a = 1) = {±2ei±e2±3e3}. (6.25) 



Figure 1 




On Figure 1 we have given the space configuration of the first particle whenever the system 
is in the state = |p; 3, 0, 8; 1, ...). The thick dots are the nests. There are 8 nests where 
the first particle can be accommodated. The coordinates of each vertex (= nest) are clear 
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from figure 1: A = (2,1,3), B = (-2,1,3), C = (-2,-1,3), D = (2,-1,3), etc. It is 
however impossible to predict which is the nest the particle is going to occupy in each 
individual experiment. 

Let Px, X = A,B, C, D, E, F, G, be the probability the first particle to be accom- 
modated in the nest X. Then proposition 6.4 asserts that with equal probability 1/2 the 
particle will be measured to have a kth coordinate ^/n^ + and with the same probability 



the kth coordinate will be — V^fc + ^4- Therefore, see figure 1, 

k=l, Pa + Pd+Pe + Ph = Pb + Pc + Pg + Pf = 1/2, (6.26a) 

k = 2, Pa + Pb + Pf + Pe = Pc + Pd + Pg + Ph = 1/2, (6.26b) 

A; = 3, Pa + Pb + Pc + Pd = Pe + Pf + Pg + Ph = 1/2, (6.26c) 

Pa + Pb+Pc + Pd + Pe + Pf + Pg + Ph ^ 1- (6.26d) 



Clearly, equations (6.26) are not enough in order to determine all probabilities Pa,---,Ph- 
There are 7 equations for 8 undeterminate. We come back to this problem in section 7, 
where the equal probability Pa = ■■■ = Ph will be proved (see propositions 7.3 - 7.5). 

The picture on Figure 1 rises many questions. The main question for us is whether 
r(|p; n), a) includes all nests where the a-th particle can be accommodated. Such a ques- 
tion make sense. The nests obtained so far are based on Eq. (6.19). Clearly, the corre- 
sponding nests can be only at the vertexes of a parallelepiped with each edge being parallel 
either to ei or to 62 or to 63. This observation suggests to search for possible nests based 
on any other triad e^, e^-, 63 of orthogonal unit vectors. Let f^;^, ^aS coordinate 
operators along e'^, 62, 63, respectively. We shall see that (r^^)^, (^^2)^? ('^as)^ commute. 
If it happens in addition that 

^ai|p; ^i, - , ••) = r2i\p; .., nj, na+3, ..) (6.27) 

then repeating the arguments from above, one would conclude that the points 

±^e;, ±y^e'2, ±^e'3, z = l,2,3, (6.28) 

are also nests for the ctth particle. This time the nests would be (measured to be) in the 
vertexes of a parallelepiped, which edges are parallel to either e'^,, 62, or 63. Therefore, 
depending on the orientation of the triad, e'^, 62, 63 the nests (6.28) could be new. We 
wish to underline that the result is certainly independent on the choice of the basis. The 
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nests (6.28) can be written in the basis (6.8), in the basis 61,62,63 or in any other 3D 
basis. From a technical point of view however the choice of the basis may be relevant. 

We find it convenient to consider all possible frames obtained by rotations from 6. 
Any such rotation is determined by a real orthogonal 3x3 matrix g = (gij), i,j = 1,2, 3, 
i.e., gg^ = 1 denotes a transposition) with determinant 1. The set of all such rotations 
constitute a group, the group SO {3). To each g put in correspondence a linear operator Tg 
rotating the basis as follows: TgCfc = {Gg)k- One verifies that the map g ^ Tg determines 
a representation of SO {3) in 3D. 



5(l)Jg(2) = Jg(l)g(2), J- E - ^ , 



(6.29) 



where E is the 3x3 unit matrix. 

Denote by 6(^) = (6(^)1, 6(^)2, ^{g)z) the orthonormed frame obtained from the initial 
one (6.8) after a rotation g\ 

3 

TgCk = e(^)fc = Gigik, k = 1,2, 3, or in matrix notation 6(^) = eg. (6.30) 

We shall parameterize the matrices g e S0{3) (and hence the frames (6.30)) with the 
Euler angles a, /3, 7 as in [56]: 



cos a cos P cos 7 


— cos a cos (3 sin 7 


cos q: sin /3 \ 


— sin a sin 7 


— sin a cos 7 




sin a cos /? cos 7 


— sin a cos P sin 7 


1 sin a sin P 


+ cos q: sin 7 | 


+ cos a cos 7 




— sin f3 cos 7 


sin P sin 7 


cos P / 



(6.31) 



The domain of definition of the Euler angles in (6.31) is 

< a < 27r, 0<p<TT, < 7 < 27r. 



(6.32) 



Different triples {a, P, 7) define different rotations apart from the cases P = when = 
a' + 7' defines one and the same rotation and P = tt when a — 7 = a' — 7' corresponds 
also to one and the same rotations around 63. 

Each rotation g{a, P, 7) can be represented as a sequence of rotations around 61, 62, 63. 
One such possibility (which reduces to rotations only around 62 and 63) reads: 



g{a, p, 7) = 5r(63, a)g{e2, P)g{e3, 7), 
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(6.33) 



where g{ek, (fi) is a rotation around on angle (p. Explicitly 

/I 

g{ei,(f)=\o cos(fi - sinip \ = e~"^^\ -si = i(e32 - 623), (6.34a) 
y sin (fi cos (f 

cos (p sin ip 

g{e2,<p)^{ 1 |=e-''^^% S2 = i(ei3-e3i), (6.346) 
- sin cp cos (p ^ 

cos 99 — sin (p \ 
5f(e3,(yi')= ( sinv? cos(p = e~''^''^ S3 = i(e2i - 612). (6.34c) 
1/ 

Here e^j are the 3x3 matrix units. Hence g{a, (3, 7) can be written also as 

g{a,P,^) = e-^"*3e-i^s2g-i7«3_ (635) 

Each state space V{N,p) is an sl{3\N) module and therefore it carries also a repre- 
sentation of the physical SO (3) group with generators (of the so(3) subalgebra) 

5'i = 1(^32 — -E23), <S'2 = i(-Ei3 — E31), ^3 = i(£^2i — -^12)- 

To each rotation g{a,(3j'y) there corresponds a "rotation" in the state space V{N,p) by 
an unitary operator U {g{a, P, 7)): 

$ ^ = J7(£?(a,/3,7))$, ^eV{N,p). (6.36) 

where in agreement with (6.35) 

U{g{a,l3,-i) = e-*«^3e-i/3S2g-i753_ (5,37) 

In particular the operators tj{g{ek, (p)), k — 1,2,3, corresponding to rotations g{ek,^) 
around on angle (p read: 

U{g{ek, ^)) = e-"^^'' , k= 1,2,3. (6.38) 

To each rotation g there corresponds also a "rotation" in the algebra of the observables, 
induced via the transformations U (g) of the state space. Indeed, consider the observable 
L and let be a normed to 1 state, which is a linear combination 



* = ociipi + a2ip2 + ••• + ocpipp, 
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of an orthonormed set of eigenstates ipi,ip2, ■ ■ ■ , ipp of L: Lipi = Xiipi. Then 

(L)* = (*,L*) = laipAi + ... + |a„pAn, 

and the superposition principle asserts that Icujp = \{^, Lipi)]'^ gives the probabihty of 
measuring the eigenvalue of the observable L. Therefore this probability cannot depend 
on the choice of the coordinate frame. Based on this Wigner proved a stronger statement 
[57]: for any two states $ and ^ the matrix element ($,L^f) should be invariant under 
any rotation g of the basis, namely 

($, L*) = {U{g)^, L{g)U{g)^). (6.39) 

Clearly this is the case only if 

L{g) = U{g)LU{g)-\ (6.40) 
For brevity we denote the above unitary transformation as V{g): 

V{g)L = ij{g)LU{g)-\ (6.41) 

Evidently 

V{gi92) = V{9i)V{92), and V{g-^) = V{g)-\ (6.42) 

Hence the map g — > V{g) defines a representation of 50(3) in the algebra of the observables 
(considered as a linear space). 

Let now the operators under consideration be the coordinate operators for particle a: 
fa = (^ai) ^a3) '^as)- How do they transform under rotations? 

Proposition 6.5. The transformation relations of the operators (fai, fa.2i Tas) under global 
rotations g are the same as for the frame vectors: 

3 

r{g)ai = U{g)raiU{g)~^ = ^rajgji, z = 1,2,3, (6.43) 

or in a matrix form 

r{g)a = {rag). (6.44) 

For a proof see Appendix A. 
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The physical interpretation of r{g)ak is the same as in the canonical QM: f{g)a,k is 
the coordinate operator of the ath particle along e{g)k- 



Since any unitary transformation preserves the commutation relations, from (6.43) 
one concludes that the operators r{g)'^j^ commute: 

[r{9)li,r{g)%]=0 y a,P=l,..,N, i,j = l,2,3. (6.45) 

Our next task is to determine the dispersion of r{g)'^ as a function of g whenever 
the system is in a fixed state i.e., Disp{r{g)'^iS)\p.n)- Then the matrices g for which 

the dispersion vanishes simultaneously for k = 1, 2, 3 will determine possible nests for the 
particle a under consideration. 

From the transformation of the basis \p;n) under the action of f{g)a,k 

f{g)c.,k\p-, ni: n2, ng; .., n«+3, ..) = (-i)-i+"2+-3+n„+3-i 

3 

XI {djkie'^'^*^ {uj + l)na+3\p; ..,%• + !,..; ..,na+3 - 1, ••) (g ^g) 

+ e~'^'^*y^ (1 - na+3)nj\p; .., nj - 1, ..; .., n^+s + 1, ••))) 

one concludes that the mean value of each coordinate operator f{g)a,k in the state \p; n) 
vanishes: 

{r{9)a,k)\p;n} = 0- (6.47) 

Again from (6.46) and the circumstance that all r{g)a,k are Hermitian operators one 
computes the mean square deviation of each coordinate operator r{g)a,k in the state |p; n): 



{r{9)l,k)\p;n) = (^(i/)a,fc|p;n),f(^)a,fe|p;n)) = X^ffe(ni + na+3), k = 1,2,3. (6.48) 
In view of (6.47) the above relation yields actually the dispersion 

3 

'Disp{r{g)a,k)\p;n) = 3iki.i^i + i^a+s), k=l,2, 3. (6.49) 



i=l 



oi f{g)a,k the state \p; n). 

Note that for a state ip = \p; n) such that ni = 77-2 = ns = na+3 = the dispersion 
(6.49) vanishes, 

Disp(f(^)a,fe)^ = 0. (6.49a) 
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Moreover such a state always exists if p < N — 3. 

We proceed to derive an analogue of the uncertainty relations for the position and the 
momentum operators. Choose an arbitrary direction in the 3D space determined by a unit 
vector n = n(6', (f) with spherical coordinates 9 and (p: 

n = h{d,(p) = (ni,n2,n3) = (sin^cos<^, sin^sin(^, cos(^). (6.50) 

Let go be the rotation matrix (6.31) with Euler angles a = (f, P = $ and an arbitrary 
1-90 = 9o{<P,d,-f). Then 

n(^,^) = e(firo)3. (6.51) 

Consequently the projection operator r{h)a = r{9, (p)^^ of the position of the ath particle 
along n coincides with the coordinate operator f{gQ)o,,z'- 

r{h)a=f{0,(p)a = f{gQ)a,3- (6.52) 

The last results together with (6.49) allows us to compute the dispersion of r{9, (fi)a for 
the cuth particle in the direction h{9, ip) whenever the system is in the state |p; n): 



D\SY>{{f{9,ip))a)\p-n) = Disp(f(5fo)a,3)|p;n) = ('^i +na+-i)cos^^ sin'^9 
+ (n2 + na+^)sin^(p sin'^9 + (ns + na+^)cos^9. 

From the inequality nj- + n^+s < p, we deduce 



(6.53) 



Disp((f(^,(^))«)|p;n) <p. (6.54) 

The last expression holds for any basis vector |p; n) from the state space W{N^p)^ for any 
direction n(^,(^) and for any particle a. Skipping these labels we can write Disp(f) < p. 
Taking into account that p^ ^ = ^ one also concludes that Disp(p) < p. Thus, for the 
standard deviations we finally obtain: 

Ar < Ap< (6.55) 

On Figure 2 we have given the standard deviation of the first particle whenever the 
oscillator is in the state \p; 3, 0, 8; 1, 77-5, njv+s)- Observe that the standard deviation is 
within a parallelepiped, the length of each edge of which is less than -\/l2 since the minimal 
possible value of p is 12. 
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Figure 2 




The conclusion is that the uncertainty of the position and of the momentum along any 
direction, for any particle and for any state is less than Consequently for the analogue 
of the uncertainty relations in QM we have 

AfAp < p. (6.56) 

Then for the initial position and momentum operators, see (6.9), eq.(6.55) yields: 

AW<^;^. (6.57) 

The above equation is very different from the corresponding uncertainty relations in 
QM, for instance AxAp^; > In particular the increasing of the accuracy of the position of 
a particle does not force the uncertainty of the momentum of the same particle to increase. 

We proceed to evaluate the dispersion of r(5f)^ ^ for any g in the state \p\ n). For this 
purpose it suffices to compute the transformation of the basis |p; n) under the action of 
f(fir)^ fc 9 k = 1,2, 3. Here is the result: 

r{9)l,k\P; = {alki^i + ria+s) + 92kin2 + n^+s) + glki^s + ^a+s)) \p; n) 

+ 9ik92k{Vi^i + l)^2|p; n)i-2 + a/ (n2 + l)ni\p; n))_i,2) 

. . (6.58) 

+ 9ik93kW{ni + l)n3|p;n)i -3 + V (na + l)ni|p; n)_i,3) 

+ 92k93k{Vi^s + ^)n2\p; n)-2,3 + \/ {n2 + 1)^3 |p; n)2-3) 
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This result is in agreement with (6.48), since {\p; n),r{g)'^ n)) gives the same expression 
for (f(^)^ A;)|p;n>- Moreover (6.58) helps to compute (r(^)^ k)\p;n) using the hermiticity of 

{■f{9)t,k)\p;n) = (b;n),f(^)^^fc|p;n)) = {r{g)l,k\p;n'),r{g)l^k\p;n')) 

(6.59) 

+ 9lk9lki'^^i^2 + ni+ 712) + gli^glkC^niUs + ni + ns) 
+ 9lk9iki'^n2n3 + n2 + ns). 

Hence for the dispersion of {r{g)a,kY = ^(fi')a k state \p; n) we finally obtain 

Disp(r(£f)^^fc)|p.„) = {r{g)i^k)\p;n) - {r{9)l,k)'\p-n) 

= fi'lA;fi'2fc(2?^l?^2 +ni+ 712) + glkgikC^nin-s + rii + m) (6.60) 
+ 5'3fe5'2fc(2n2n3 + n2 + n3), /c = l,2,3, Q; = l,...,iV. 

which can be written also as 

3 

Bisp{r{g)l^k)\p;n) = ^ aikQ^kC^n^nj + Ui + Uj) . (6.61) 
i<j=i 

Then the expression for the standard deviation of f{g)^ ^ in the state |p;n) read: 

^{r{9)a,k)\p;n) = (i/lfei/2fc(2nin2 + Tli + n2) 

.1/2 (6-62) 
+ gikglkC^nms + ni + ns) + glj,glk{2n2n3 + n2 + ns) j . 

The problem to solve now is to find all different reference frames e(^) for which 
the dispersion (6.61) vanishes simultaneously for all three values of A; = 1,2,3. If e(^) = 
(e(^)i, 6(^)2, e(^)s) is one such frame, then according to Conclusion 6.2 the state |p; n) will 
be a common eigenstate of r{g)a,i: ^{9)0,2 ^^^l r(^)^ 3. This means that the nondiagonal 
terms in (6.58) have to vanish so that: 

^(^)a,feb; = ((^ifcK + ^a+3) + 92kin2 + ria+s) + ^|fc(^3 + ria+s)) \p; u) , (6.63) 
The eigenvalues 



ri9)l,k = fifeK + ^a+s) + 92ki'r^2 + ria+s) + glkins + Ua+s), 
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(6.64) 



of 'r{g)'^ k = 1, 2,3, are the squares of the admissible coordinates of the ath particle in 
the frame e(^). Let us summarize. 



Conclusion 6.2. Let the system be in the state \p',n). If the dispersion Disp(f(^)^ k)\p;n) 
vanishes for a certain g and for all k = 1,2, 3, i.e., 



(6.65) 



^^Mr{9)l,k)\p;n) = 9fk9iki'2nin2 + ni + n2) 
+ i/ifef3fe(2^i^3 + ni + ns) + glkdlki^'T^^ns + + ns) = 0, 

then \p;n) is an eigenstate of r{g)'^ ^, 

r{9)l,k\p-, n) = r{g)l^k\v^ n), (6.66) 

with eigenvalues 

^(i/)a,fe=^lfe(^l+^a+3)+^ife(^2+na+3)+^L(n3 + na+3), k = 1,2,3. (6.67) 

In such a case 

T{\p; n)),a, g) = {r{g)oc,M9)i + r{g)a,2^{g)2 + r(^)a,3e(^)3} (6.68) 

determines admissible places, i.e., nests for the ath particle. 

The validity of the above proposition can be verified also directly. To this end write 
the dispersion (6.60) as follows: 

D{f{g)l^k)\p;n} =9lk9lk{ni + 1)^2 + ^L^ifc(n2 + l)ni 

+£'ife£'3fe(^i + 1)^3 + glkglkin-s + 1)^1 (6.69) 

+f2fcf3fc(^2 + 1)^3 + 9lk9lkin^ + l)n2. 

Since the dispersion (6.69) is a sum of nonnegative terms, it vanishes only if every term 
vanishes, i.e. if for any i < j = 1,2,3, 

9ik9jk{ni + l)nj = and g1k9%{jij + l)ni = 0. 

But then 



9ik9jk\J{ni + l)nj = and gikgjk\J {nj + l)ni = 0, i < j = 1,2,3 

and therefore all off diagonal terms in (6.58) vanish. Hence (6.66) and (6.67) hold. 
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The conclusion 5.7 is clear. What is not so clear is how many if any are the new 
nests in addition to (6.24). We shall answer this question first on the example, namely 
one-particle 3D oscillator with p= 1. 



Example 6.2. Consider one-particle Wigner oscillator in a representation p = 1. The 
state space V{N = l,p = 1) is 4-dimensional with a basis 

^1 = |p = 1; 1, 0, 0, 0), ^2 = |p = 1; 0, 1, 0, 0), = b = i; o, o, i, o), ^4 = b = i; 0, 0, 0, 1). 

(6.70) 

la. Take first = \p = 1;0, 0, 1;0). The requirement the dispersion Disp{f{g)1)^^, 
g — g{a, (3, 7), to vanish reads, see (6.60): 

Disp(r(fi')fc)v'3=^u5'ifc + 5'2'fc5'ifc = 0, A; =1,2,3. (6.71) 

Since the system consists of only one particle, we have suppressed in (6.71) the index a. 

The simplest solution of (6.71) corresponds to a = /? = 7 = 0, namely to g being a 
3x3 unit matrix, g = 1. In this case 

fl(p3 = 0, fl(p3 = 0, fjifis = (p3. (6.72) 

Then according to Conclusion 5.7 

r{\p = 1; 0, 0, 1; 0), a = 1, ^ = 1) = {±63} = (0, 0, ±1), (6.73) 

i.e., there are two nests A with coordinates (0, 0, 1) and B with coordinates (0,0,-1) where 
the particle can be registered to be (see Figure 3). 

Apart from g = 1, the Eqs. (6.71) have also other solutions. It takes some time to 
find all of them (see corollary B.l for more details): 

a) e {0, tt}, a, 7 —arbitrary, (6.74a) 

b) /? = 7r/2, 7 e {0, 7r/2, TT, 37r/2}, a - arbitrary, (6.746) 

Do they lead to new nests? No, they do not. It turns out that any of the above solutions 
corresponds to the same picture shown on Figure 1. Let us illustrate this on an example 
with /? = 7r/2, 7 = and a arbitrary. In this case 

(0, —sin a, coso: 

0, cos a, sincK | , (6.74) 

1, 
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Figure 3 




and (6.58) yields: 

r{g)lv^3 = ^3, r{g)lip3 = 0, f{g)lip3 = 0. (6.76) 

Therefore 

r{g)j = 1, r{g)l = 0, rig)l = 0, (6.77) 
so for the nests we obtain applying (6.68) 

T{\p = 1; 0, 0, 1; 0), a = 1, ^73) = {±e(^7)i}. (6.78) 

But e{g)i = J2j ^j9ji = ^s? i-^-, we obtain the same nests as in Figure 1: 

T{\p = 1; 0, 0, 1; 0), a = 1, <7 = 1) = r{\p = 1; 0, 0, 1; 0), a = 1, gs) = {±63}, (6.79) 

The proof for the other cases in (6.74) is similar and it gives all of the time the nests 
as shown on Figure 1. 

lb. In an analogues way one finds that the space configuration of the state \p = 
1; 1,0, 0,0) (resp. \p = 1; 0,1,0,0),) corresponds to a picture with two nests ±ei (resp. 
±62). 

Ic. The above results suggest that the nests r((^fc), corresponding to (7 = 1, see 
(6.73), contain already all nests where the particle can be accommodated. Is this the case? 
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No, it is not. In order to show this we consider the last basis state (^4, corresponding to 
77-1 = 77-2 = ns = and 71-4 = 1. 

Evidently, see (6.60), the dispersion Disp(f(^)^ jj^^ = for any g. Then, as it should 
be, <^4 is an eigenstate of 'r{g)1 for any g. Indeed, setting in (6.58) ni = n2 = ns = 0, one 
has 

r{9)l^4 = {gfk + glk + 9lk)^4- (6.80) 

The matrix g is an orthogonal matrix, gg^ = 1, and therefore gfj^ + g2^. + g^]^ = 1. Thus, 
for any g, 

fig)l(p4 = (P4, /c=l,2,3. 

Since f'{g)1 is square coordinate operator along e{g)k and the basis e(g) is orthonormed, 
the conclusion is that any point on a sphere with radius V3 is a nest. 

We see that the properties of the state (f^ are very different from the properties of the 
other three basis vectors. In particular if the system is in this state, the particle can be 
observed in every point of the sphere with radius \/3. 

Passing to the general case we divide all states into three nonintersecting classes. 

Class I. All states \p; n) with ni — n2 — — 0; 

Class II. All states \p; n) for which two of the integers tt-i, ri2, do not vanish. 

Class III. All states \p;n) for which two and only two of the integers ni, 77-2,77.3 vanish. 

1. Properties of the states from Class I 

Observe first of all that the Class I is not empty only ifp<N. The latter stems from the 
observation that the sum of the fermionic coordinates Uf of any state \p; n) cannot exceed 
AT, whereas rii + ... + un+s = p. Therefore if p > N, then at least one of the bosonic 
coordinates ni,n2,n3 of \p;n) cannot vanish. 

la. The next property is almost evident from the results just proved (see part Ic in 
the example 6.2). 

Corollary 6.3. If the system is in a state \p;0, 0,0; ...,1^+3, ■■) from Class I, then any 
point on a sphere with radius y/3 is a nest for athe particle. 
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lb. The state = \p; 0, 0, 0; tt-q+s = 0, ..) from Class I is of particular interest. One 
verifies that 

Hai^ = Sc^k^ = Rak^ = Pak^ = 0, k = 1,2, 3. (6.81) 

Hence, 

Corollary 6.4. A state 

\p; 0, 0, 0; n4, na+2, Oa+3, ?^a+4, nN+s) (6.82) 

corresponds to a space configuration of the system when the ath particle "condensates" 
onto the origin of the oscillating system with zero energy, zero momentum and zero angular 
momentum. 

The property that some of the oscillating particles can condensate onto the origin 
with zero energy exhibits another difference with the conventional case: in canonical QM 
the ground energy of any 3-D free harmonic oscillator is never zero, it cannot be less ^u>h. 

2. Properties of the states from Class II 

We have already indicated in proposition 6.4 that if the system is in the state \p, n) 
then the collection of points 

r(|p; n))« = {±y/ni + na+3 ei ± A/nT+rWhS &2 ± a/^s+^Ws &^}\ (6.83) 
are nests for ath particle. Now we prove a stronger statement. 

Proposition 6.6. If the system is in a basis state \p; n) from Class II, then the nests 
(6. 88) are the only nests for the ath particle. 

For the proof of this relatively long proposition see appendix B. 

3. Properties of the states from Class III 

It remains to investigate the space structure of the basis states from the Class III, 
namely all those basis states for which two and only two of the bosonic coordinates 
ni,n2,n3 of \p;n) vanish. 
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Let us consider for definiteness a state \p; 0, 0, 77-3, .., n^+s, ...). In this case, see (6.60), 
the condition the dispersion of to vanish reduces to 

{9lk9ik+92k9ik)n3 = 0, (6.84) 

and since ris ^ 0, (6.84) is equivalent to 

9lk9lk = 0, 9lk9lk = ^- (6.85) 

We have aheady found all solutions of Eqs. (6.85), see corollary B.l (in appendix B). If 
the (7— matrix is one such solution, then the state \p\ 0, 0, ns, na+s, ..) is an eigenvector 

of ^(^)afc' 

^{9)lck\P'i 0> 0> ^3, na+3, ••) = r{g)'ij,\p; 0, 0, ns, n^+s, ••) (6.86) 

with 

^{9)ak = alk^is + na+3- (6.87) 

Assume the system is in the state \p; 0, 0, 72,3, n^+s, ..) from Class III. Also in this 
case the results depend essentially on the value of Uct+s- 

3a. Let n^+a = 0. Then the ath particle has only two nests, namely ±^^63. 
The proof of this result is essentially the same as in example 6.2, part Ic, so we omit 
it. In a similar way one establishes the space structure of the ath particle in the states 
|p; 0, 77-2, 0, Oq,+3, ..) and |p; ni, 0, 0, Oq,+3, ..). The results are collected in the next 
proposition. 

Proposition 6.7. The correspondence state - space structure for the ath particle read: 

|p;0,0,n3,...,0Q,+3, ..) ^ ±v^e3, ns ^ 0, (6.88a) 
|p;0,n2,0, ...,0q+3, ..) ^ ±^/n^e2, n2 7^ 0, (6.886) 
|p;ni,0,0, ...,0o,+3, ..) ^ ±y/rHei, m 7^ 0, (6.88c) 

3b. The case with Ua+s = 1 is more involved. Here is the result. For the proof see 
Appendix C. 

Proposition 6.8 

1. Whenever the system is in a state \p; ni, 0, 0, .., 1q+3, ..) all nests of the ath particle are 
r(|p;ni,0,0, ...,la+3, ..), a) = {CiVm + 1 ei + (cos ck - sin Q;)e2 

(6.89) 

+ (sina + cosQ!)e3 I en e R, ^1 = ±1, }, ni 7^ 0. 
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2. All nests of the ath particle whenever the system is in a state \p; 0, n2, 0, .., la+3, ••) ore 

r(|p;0,n2, 0, la+3, ■■), a) = {(cos a - sinQ;)ei + C2\/n2 + 1 e2 

(6.90) 

+ (sin a + cos Q;)e3 I en e R, ^2 = ±1, }, n2 7^ 0. 

3. All nests of the ath particle whenever the system is in a state \p; 0, 0, n^, ■■, la+3: •■) (^^^^ 

r(|p; 0,0, 77,3, ...,1q+3,..), a) = {(cosQ! - sinQ!)ei 

' : (6.91) 

+ (sin a + cos a)e2 + ^sVns + 1 63 | a e R, ^3 = ±1, }, ^3 7^ 0. 

From the results obtained so far we can draw already some conclusions about the 
collective properties of the system. First we observe that the dispersion (6.60) for the ath 
particle does not depend on a. This leads to the following conclusion: 

Corollary 6.5. Let the system be in an arbitrary basis state Then if the dispersion 

Disp(f (g)^^^ j^)|p.n), k = 1,2,3, see (6.60), vanishes for one particular particle a^, then it 
vanishes for all particles. Therefore if\p;n) is an eigenstate of f{g)\^^^^, ^(fi')ao,25 ^(fi')ao,3 
for one particular particle a^, then it is an eigenstate of r{g)'^ i, r{g)'^2, '^{9)a 3 f^"^ 
all particles a = 1,2,...,N. Consequently, all observables r{g)ai, 'f'{9)a,2} '^i9)a,3' ^ — 
1,2,...,N can be measured simultaneously whenever the system is in the state \p',n). 

Observe next that whenever the system is in a state |p; n) the coordinates of the nests of 
a particle do not depend explicitly on a. The dependence is indirect, via na+3, which 
can take only two values, or 1. As a consequence one concludes: 

Corollary 6.6. Given a basic state \p',n). All particles ai, 0:2, o;/- for which Ua^+s = 
^^02+3 = ••• = iT'ak+3 = have common nests; the nests of the rest of the particles, namely 
those for which noi^^^ = ^01^+2 = ••• = = 1 ^^•s^ coincide. 

Consider for instance a 6— particle oscillator in a state (pi = |p = 4; 0, 0, 1; 0, 0, 0, 1, 1, 1) 
(from the Class III). Then according to proposition 6.7 the first three particles #1,2,3 
have two common nests: 63 and —63. The space configuration of particles 4,5,6, is very 
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different. Similar as on Figure 3, each particle is accommodated somewhere on two circles 
with radius -\/2 around 2;— axes, which are on a distance -\/2 above or below the xOy. 



7. Physical properties - angular momentum, parity and probability distribu- 
tions 

Here we discuss as a first step the properties of the angular momentum of the oscillator 
system. The results are very different from the angular momentum properties of both the 
conventional oscillators and the WQSs studied so far. The most unusual new feature is that 
the operators of the projections of the angular momentum are the same for all particles, 
see (3.28), and they coincide with the generators of the algebra so(3) of the rotation group. 
As a second step we introduce another important physical observable, namely the parity 
P of the states and finally we use it in order to show that any particle occupies with equal 
probability anyone of its nests. 

Explicitly the physical so(3) generators, defined as operators in W{3\N) read: 

Si=i{b+b^ -b+b^), ^2= 1(6+63- -6+6r), ^3 = i(62+6r-6+62-). (7.1) 

Then the angular momentum projections of the ath particle are 

Mai^ji^S,, ^=1,2,3, (7.2) 

whereas for the components of the angular momentum of the entire oscillator one has 

M,. = -^5,-. (7.3) 

As a result the oscillating particles behave as if they were charged particles in a strong 
magnetic field: the angular momentums of all particles are parallel to each other. 

The operators (7.1) are not diagonal in the basis (4.2). We proceed to introduce a 
new basis which diagonalizes ^'3. To this end consider the unitary matrix 

/7I 71 0\ 



7i 71 
10 



(7.4) 



\ Iat/ 

where Iat is an iV— dimensional unit matrix. Let c+ = (c+, c^, . . . , c^^^). Denote by 
the Hermitian conjugate to G matrix. Then the operators 

c(G)+ = (c+G),, c{G)- = {G+c-)u z = 1,2,.., AT + 3, (7.5) 
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satisfy the conditions listed in (4.1). More explicitly, 

1. ciG)t ^ B+ = ±{b+-ibt), c{G)t ^ B+ = ±{bt + ibt), c{G)t ^ B+ = 6+ (7.6) 

c(G)r = Sr + i^>2 ), c{G)+ ^ B- = ^{K- '^K), c{G)^ ^ ^a" = ^3 , (7-7) 

are Bose operators and odd elements (as linear combination of odd elements), 

2. c(G)J = F±=/±, A = A,5,...,N+3, (7.8) 

are Fermi operators and even elements. 

3. The Bose operators anticommute with Fermi operators. 

An immediate consequence of the above results is that all vectors 

, , {c{G)tr^{c{G)tr^ . . . (c(G')+^3)--+3 
|p;ni,n2,...,njv+3) = / , , , 10), 



{B+r^{B+r^{btr^{f+r^ . . . (/+)"-+ 



■|0), (7.9) 



^/nl\n2\n3\ 
where 

ni,n2,n3eZ+, n4, ns, . . . , tiat+s e {0, 1}, a,nd ni + ... + un+s ^ P, (7-10) 

constitute an orthonormed basis in 1^(3|A^). The transformation of the new basis under 
the action of the CAOs Bf, S^, B^ and the Fermi CAOs is the same as in (4.4) with bf 
replaced by B^, i.e. 



B+\..,ni,..) = Vni + l\..,ni + l,..), i = l,2,3; (7.11a) 

B~\..,ni,..) = ^/r^i\..,ni-l,..), i = 1,2,3; (7.116) 

f+\.., m, ..) = (-l)-i+-+"^-ix/r^|.., m + 1, ..), i = 4, 5, .., iV + 3; (7.11c) 

f-\..,ni, ..) = (-l)-^+--+"*-i V^|..,ni - 1, ..), z = 4, 5, ..,iV+ 3; (7.11d) 
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Taking into account that 



bt = ^{Bt + B+\ bt = ^{Bt-B+), bt = B+, (7.12a) 

b- = l={B-+B-), b^ = ^{B--B-), b^=B-, (7.126) 
we obtain from (7.1) 

^1 = Kbtb2 - btb^) = -^WBT - BtB^ + B+B^ - B+B^), (7.13a) 

^2 = Kbtbs - b+b^) = -^i-B^B^ - B+B^ + B+B^ + B+B^). (7.136) 

^3 = Kbtbl - b+b^) = B+B^ - B+B'. (7.13c) 

Then 

S+ = V2{B+B^ - B+B^), = V2{B+B^ - B+B^), (7.14) 

and 



S+\p;ni,n2,n3, . . . ,njv+3) = •\/2(n3 + l)ni|p;ni - 1,77,2,77,3 + 1, . . . , tijv+s) 



- ^/2{n2 + l)n3\p; ?7i, ?72 + 1, ^3 - 1, . . . , ?7Ar+3), (7.15a) 

S'-Ip; 771,772, 773, • • • , ^^-AT+s) = -\/2(n7TI)n^b; Tli + 1,712,713 - 1, . . . , 71JV+3) 

- y/2{n3 + l)n2\p; ni,n2 - l,7i3 + 1, . . .,7iiv+3), (7.156) 
The operator ^3 is diagonal in the basis |7ii, 712, • • • , 7iAr+3): 

5'3|p;71i,7l2,...,7liv+3) = (712 - 71i) |p; Tli , 7l2, . . . , 7lAr+3) • (7.16) 

Hence Ss = i(6^6^ ~^i^^2 ) diagonal in the basis (7.9), but written in terms of the initial 
CAOs (7.12) namely 

(-^(^>^-i^>^))"H-)5(^>^+i^>^))"^(^>^)"^(//)"^•■■(/i^+3)""+^ , 

|p;71i,7l2,...,7lAr+3) = , , , , 10) 

Vni!n2'7i3' 

(7.17) 

We call the basis (7.9) (resp. 6.17) S^-basis and abbreviate 



|p;7ii, ...,7iiv+3) = \p;n). 
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(7.18) 



For the so(3) Casimir operator 



s^ = sl + si + si = s+s_ + si- Ss (7. 19) 

we find 

S^|p;ni, ..,nN+3) =(2(ni + l)?i3 + 2(?i3 + l)n2 + (n2 - ni)^ - n2 + ni)|p; ni, ..,nN+3) 



-2-y/ (ni + l)(n2 + 1)^3(713 - l)|p;ni + l,n2 + l,n3 - 2,n4, ..,niv+3) 

2-\/ (^3 + l)(n3 + 2)77,177,2 |p; 7li - 1, 7l2 - 1, ^3 + 2, 7l4, .., 7ljv+3) 

This result is not unexpected. The so(3) Casimir operator is not proportional to the unity 
because the corresponding representation is not irreducible (and is not a direct sum of 
irreps with the same signature). For further use we formulate a 

Proposition 7.1. All vectors \p',ni,n2,n3,n4, ...jUn+s) with one and the same nt, = 
ni + 7i2 + 7i3 and fixed 714, un+s are vectors from V{N,p, rib, ^4, un+s) and therefore 
they have one and the same energy (5.4). 



Indeed, 



y/2 ^ \ y/2 

ni 712 

where c{ni,n2, k, q) are numbers. Therefore, 

|p;7ii,7i2,...,niv+3) = ; , , , 10) 

ni n2 

= ^^d{ni,n2,k,q)\p;ni + n2 -k-q,k + q,ns, ...un+s), 

k=0 q=0 



(7.20) 



where again d{ni,n2,k,q) are numbers. Clearly all vectors in the RHS have one and the 
same rib and therefore they represent states with one and the same energy. 
Next we simplify some of the notation: 

rib = b = rii + n2 + ns - the number of the bosons in a state \p; n) , 
7i/ = / = 7i4 + ... + 71JV+3 - the number of the fermions in \p; n). 
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We are now ready to describe the so(3) structure of V{N,p). The first step of the 
problem was aheady carried out (see (5.9)): 

p 

V{N,p)= Vi{N,p,nb)<»V2{N,p,nf=p-nb), (7.21a) 

rib =max(0 ,p—N) 

where Vi{N,p,nb) is an irreducible ^/(3) module, see (5.12), and V2{N,p,nf), see (5.14) 
is an irreducible gl{N) module. Another way to write (7.21a) is 

p 

V{N,p) = 0(iV -p + nb)Vi{N,p,nb) (8) V2{N,p,nf =p- m,), (7.216) 

where Q{x) = for a; < and 1 for a; > 0. 

Since so(3) is a subalgebra of gl{3), the rotation algebra transforms only the bosonic 
part Vi{N,p,nb) of Vi{N,p,ni)) (8) V2{N,p,nf = p — Uh). Hence in order to determine the 
angular momentum structure of the system in V{N,p) one has to decompose each gl{3) 
module Vi {N, p, b) along the chain 

gl{3) D so{3) D so{l). (7.22) 

For the ladder representations of gl{S), which we consider, the problem was solved in [58] 
directly for the quantum case. We shall use the results from [58], but without deformations. 
One possible orthonormed basis in Vi{N,p,b), consistent with (7.20) is 

\p;ni,n2,ns)= ^ , ? \ ^ T^^\Q), K\^) = ^^ + + = b = p - f. (7.23) 

In this basis the so(l) generator 5*3, see (7.16) p. 37, is already diagonal. The basis (7.23) 
is not however reduced with respect to so(3). 

The decomposition of Vi{N,p, b) into irreducible so(3)-modules reads [58]: 

V,{N,p, b) = ^V^{N,p,b,S), 5 = 6, 6 - 2, l(or 0), (7.24) 
s 

where Vi{N,p,b, S) is an irreducible so(3) module with angular momentum S. As an 
appropriate orthonormed basis in Vi (N, p, b, S) one can take 



v{p,b,S,Ss) = 



{b + s)\\{b - sy.lis + Ss)Ks - S3)\{2S + 1) 



X 



{b + S+l)\ 

L(S+S3)/2J (6-S)/2 ^(^3 + b-2x- 2y)l{2x + 2y)!!(2x + 2y- 2Ss)l 



...ax(o,53) .^0 (2x)n(2y)!!(2x - 2S,y.KS + - 2x)Kb - S - 2,)!! 



X \p; X + y - Ss,x + y,b + Ss - 2x - 2y), (7.25) 
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where ^'3 is the projection of the angular momentum along the 2;— axes. 

Observe that the coefficients in the RHS of (7.25) do not depend on p. Therefore if 
the relation 

6, 5, 53) = c{b,S,S3,ni,n2,ns)\p;ni,n2,n3) (7.26) 



ni,n2,n3 



holds for a certain p, then it holds for any p. 

As an orthonormed basis in V{N,p^ we take 

v{p,h,S,S3)®{j^T\hTK..{j^^3Y-^-\^), (7.27) 

where (/4)"''(/5)"'®...(/Ar+3)"'^+^ |0) with 714 + ... + nAr_|_3 — nf—p — rib is an orthonormed 
basis in V2{N,p,nf). Instead of (7.27)we shall also write 

^;(p,6,5,53)®(/4)"^(/5)"^..(/^+3)""+lO) = \\N,p,b,S,S3,n^,n5,...,nN+3))- (7.28) 

The decomposition of any basis state (7.28) in terms of the S'3— basis (7.18) follows 
from (7.25): 



(6 + S)\\{b - S)\\{S + Ss)\iS - Ssy.{2S + 1) 



\\N,p,b,S, 83,71^,715,..., un+s)) = 

L(S+S3)/2J {b-S)/2 

\^^Uo,Ss) U ' (2.)!!(2y)!!(2.-253)!!(5 + 53-2.)!(6-5-2,)!! 



(6 + 5 + 1)! 

yo-^. V(g3 + b-2x-2y)!(2x + 2y)!!(2x + 2^- 2^3)!! 

.i?o,53) h (2x)!!(2y)!!(2a;-253)!!(-S + 53-2a;)!(6-5- 

X |p; a; + y - S'3, a; + y, 6 + S'3 - 2a; - 2y; 77,4, n^, njv+3), (7.29) 



We call the basis (7.28) SO (3) -reduced basis or an angular momentum basis. By 
construction each basis vector \\N,p,b, S, S3,n4,n5, ...,nN+3)) is an eigenvector of H, 
and S3: 

H\\N,p, rib, S, S3, n4, njv+3)) = (^'^^ + ^nf)\\N,p, rib, S, S3, 77,4, njv+3)), 

S'^\\N,p,nb, S, S3,n4, ...,nN+3)) = S{S + 1)| |iV,p, n^, 5, ^3, 77,4, njv+3)), 

S'3||Ar,p,?7,5,S',5'3,?7,4, ...,nN+3)) = SsWN, p, Ub, S, S3, n4, Un+s)) ■ 

The admissible values of 77-4, ...,njv+3 distinguish between the basis states with the same 
energy, angular momentum and its third projection. Let us summarize. 
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Corollary 7.1: A reduced basis vector \\N,p,b, S, 83,714,715, ...jTIn+s)) corresponds to 
a state of the system with energy E = u!h{3p + Nb — 3b) /\N — 3\, angular momen- 
tum S, its third projection Ss and fermionic numbers n^, ...,nN+3- All different states 
\\N,p,b, S, 83,714,715, ...,nN+3)) , namely those with 

(a) b = max(0,p — A^), max(0,p — A^) + 1, ...,p — l,p, 

(b) S = b,b-2,...,l{or 0), 

(c) 83 = -S,-S + l,...,S, 

(d) fermionic numbers 714,715, ...,njv+3 such that n4 + ... + njv+s = P — b, 

constitute an orthonormed basis in the state space V{N,p). 
Let us give an example. 

Example 7.1. Let N = 2 and p = 1. The state space is 5 dimensional. The angular 
momentum basis read: 

\\N = 2,p=l,b=l,S = 1,83 = 1,0,-1,714 = 0,715^ 0)), (7.30a) 
\\N ^2,p^l,b ^0,8 = 0,83 = 0,714 = 1,715 = 0)). (7.306) 
\\N = 2,p=l,b = 0,8 = 0,83 = 0,714 = 0,715 = 1)). (7.30c) 

In terms of the S'3— basis (7.17) and the initial basis the above states read (we skip the 
common for all states labels N = 2 and p = 1): 

||6= 1,5 = 1,53 = l,n4 = 0,715 = 0)) = -|0, 1,0, 0,0) 

= --^ |1, 0, 0, 0, 0) - -^|0, 1, 0, 0, 0), (7.31a) 



= --^ |1, 0, 0, 0, 0) - -j=\0, 1, 0, 0, 0), (7.31a) 

\\b=l,8= 1,83 = 0,n4 = 0,n5 = 0)) = |0, 0,1,0,0) = |0, 0,1, 0,0), (7.316) 
\\b= 1,8 =1,83 = -1, 714 = 0)) = b = 1; 1, 0, 0, 0, 0) 

= (|1, 0, 0, 0, 0) - i|0, 1, 0, 0, 0)), (7.31c) 

||6 = 0,;S = 0,^3 = 0,n4 = 0,n5 = 1)) = |0, 0,0, 0,1) = |0, 0,0, 0,1), {7.31d) 
||6 = 0,<5 = 0,^3 = 0,n4 = l,n5 = 0)) = |0, 0,0, 1,0) = |0, 0,0, 1,0). (7.31e) 

There are 5 states. The first three of them correspond to orbital momentum 1. The 
last two states have orbital momentum 0. For an example corresponding to any N and 
p = 1,2,3 see appendix D. 

We recall, see (3.29), that the angular momentum of the particles is measured in units 
whereas for the entire system it is Mj = -^^Sj. 
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The last physical observable which we are going to consider is the parity operator P, 
called also space inversion operator. In the 3D space this operator transforms the frame 
vectors into their mirror images: 



Pefe = -efc, A;= 1,2,3. (7.32) 

In a matrix form 




(7.33) 



The matrix P is an orthogonal matrix and therefore P e 0(3). In fact [57] 

0(3) = {g, -g = gP\g G SO{3)}. (7.34) 

Our problem is to find out how P acts in the state space V{N,p). We shall use the 
circumstance that P is an element also from the group GL(3), P e GL{3). This is evident 
since GL{3) is the collection of all 3 x 3 matrices with determinant different from zero. As 
a second step we use the exponential map: [59] if a; e 91(3), then expa; e GL{3) in order 
to find (first in the defining 3x3 representation) which is the element x from the algebra 
gl{n) for which expa; = P e GL{3). 

To begin with take the element -1(^(^11+^22 + ^33) from gl{3). Then (10302) g{ip) = 
exp(-i(/?(Eii + E22 + ^33)) e GL{3). Therefore 

g{ip) = exp{-iif{Ei^ + E22 + Es3)) = ^ ^-^(£^11 + E22 + £^33)^ (7.35) 

But in the defining (3 x 3) representation En + E22 + £'33 = E is the 3x3 unit matrix E 
and E^ = E. Consequently 

00 / ■ \ fc 

9^^) = E ^-^^ = ^exp(-i<^), (7.36) 

fc=0 

and for = TT (7.36) yields: ^(tt) = £'exp(-i7r) = Ecostt = —E = P. Thus, 

P = exp(-i7r(£;ii + £;22 + ^33))- (7.37) 

The relevance of the last relation stems from the observation that it holds in any represen- 
tation of the Lie algebra gl{3) and we know the realization of the 5'/ (3) generators Ekk in 
any state space V{N,p). Indeed according to (4.7) Eij — bfbj , i,i — 1, 2, 3, and therefore 

P = exp(-z7r(6^6]- + 6^63 + htb^)). (7.38) 
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Taking into account that 6^6^ are number operators, see (4.5c), we find: 

P|p,n) = exp(-i7r(ni + n2 + n3))|p,n) = {-iy'+'^^+''^\p,n). (7.39) 

We see that the basis vectors \p; n) are eigenvectors of the parity operator and since ni + 
n2 + ns — b we conclude: 

P\p,n) = (-l)"^+"^+"3|p,n) = {-lf\p,n). (7.40) 

A number of consequences follow from (7.40). 

Corollary 7.2. Any state \p; n) is invariant under the action of the parity operator. 

Obviously = 1 and therefore P — P~^. Moreover for any two basis states \p;n) 
and \p; n') 

(Pb;n),b;n')) = (b; n), P^; n')), {P\p;n) ,V\p;n')) = {\p;n),\p;n')). (7.41) 
Corollary 7.3. P is an unitary Hermitian operator. 

The observation that the energy of a state b;n) from V{N,p) is in one to one correspon- 
dence with 6, see (5.4), yields: 

Corollary 7.4. All states \p]n) from V{N,p) which have one and the same energy have 
also one and the same parity (—1)^. In particular all states from V{N,p, b, f) have parity 
(— l)**. Consequently the parity of the state \ \N,p,b, S, 83,714,715, ...jUn+s)) , see (7.28), is 
also 

Corollary 7.5. It is straightforward to verify that 

PHP-^ = H, ^ [F,H] = 0, (7.42a) 

PrafcP"' = -rak, k = l,2, 3, (7.426) 

FpakP-^ = -Pak, k = l,2, 3, (7.42c) 

P^fcP-i = Sk, k= 1, 2, 3, (7.42d) 

Thus the parity operator P is an integral of motion, the Hamiltonian H is a proper scalar 
operator, the position and the momentum operators are proper vector operators, whereas 
the angular momentum is pseudovector operator. 
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So far we have clarified what are the possible nests for any particle whenever the 
system is in a basis state \p; n) . What we have not clarified yet is what is the probability 
the particle to occupy one of the nests. This issue will be the topic of the next discussion. 
We will show that with one and the same probability the particle under consideration can 
be in anyone of its nests. 

We begin with the states from the Class II. If the system is in a state \p; n) from this 
class, then the measurements of the coordinates of the a— th particle yield that its nests 
are on a sphere which form the vertices of a rectangular parallelepiped 

Tai = ±\/ni +na+3, = ±^n2 + na+3, = ±^n3 + na+3k = 1, 2, 3. (7.43) 
First we derive a few preliminary results. 

Proposition 7.2. Any state \p; n) from Class II is invariant under rotation on angle tt 
around x—, y— or z — axes. 

Proof. Consider a rotation of the system on an angle tt around 2;— axes. Since any state 
|p; n) is defined up to a multiplicative constant, we have to show that \p; n) is an eigenstate 
of the rotation operator exp(i7rS'3). To this end we expand the state \p;n) in the basis 
77.)}, namely via the eigenvectors of ^3. Using relations (7.12) we calculate: 

ni n2 

\p; ni, n2, ...) = X] c(ni, n2, k, q)\p; ni + n2 - k - q,k + q, ns, nN+3), (7.44) 

k=0 q=0 

where c(ni, 712, k, q) are numbers. Then from (7.16) and taking into account that exp(i7r) = 
— 1), we calculate: 

exp(i7rS'3)|p;ni,n2, ...) 

ni n2 

c{ni,n2, k, q) exp(i7r(2A; + 2q-ni- n2))\p; rii + 712 - k - q,k + q, 77,3, njv+3) . 

k=0 q=0 

= (-l)"^+-^|p;ni,n2,...) 

(7.45) 

In a similar way, a replacement in (7.12) 1^2, 2^3, 3^1, (resp 1 ^ 3, 2 ^ 1, 
3 — > 2) diagonalizes Si (resp. 82). Then repeating the above arguments one proves that 
any state \p; n) is invariant under rotation around x— or around y— axes on angle tt. This 
completes the proof. 
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Proposition 7.3. Let the system be in a state \p;n) from the Class II. Then the a—th 
particle will occupy with equal probability any one of the corresponding nests. 

Proof. Assume for definiteness that ni ^ and n2 7^ 0. There are two cases to be 
considered. 

Figures 4 




H 



a) 713 + 71q_|_3 7^ 0. The configuration of the nests is shown on the LHS of Figure 4. 
We have denoted by A the nest with all positive coordinates: 

A = + na+3, y^n2 + na+s, ^/n3 + n^+s)- (7.46) 

The coordinates of the rest 7 nests are also clear. They differ only by signs from those of 
A. 

Perform now a rotation on angle tt around 2;— axes. Clearly it will bring the system in a 
configuration shown on the RHS of Figure 4. But according to proposition 7.2 both space 
configurations shown on Figure 4 represent one and the same state. This in particular 
means that 



PiA) = P(C), P{B) = P{D), P{E) = P{G), P{F) = P{H), (7.47) 
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where P{X) is the probabihty the (a-th) particle to be in the nest X, X = A, B, ...H. 

In a similar way, performing a rotation on angle tt around axes and using again 
proposition 7.2, one concludes that 

P{A) = P{H), P{D) = P{E), P{B) = P{G), P{C) = P{F) (7.48) 

which together with (7.47) yields: 

P{A) = P{C) = P{F) = P{H), P{B) = P{D) = P{E) = P{G). (7.49) 

A rotation around axes on angle tt does not give anything new. As a next step we use the 
circumstance that the state \p\ n) is invariant also with respect to parity transformation, 
see (7.40), which yields: 

P{A) = P(G), P{B) = P{H), P{C) = P{E), P{D) = P{F). (7.50) 

The latter together with (7.48) proves proposition 7.3: 

P{A) = P{B) = P{C) = P{D) = P{E) = P{F) = P{G) = P{H). (7.51) 

b) ns + ria+z = 0. In this case the nests are only 4 and all of them are in the xOy plane. 
The proof of the equal probability given above works perfectly well also in this case and 
is even simpler. 

Let us turn to Class III configurations. 

Proposition 7.4. Let the system be in a state \p;n) from the Class III. Then the a—th 
particle will occupy with equal probability any one of the allowed nests. 

Proof. Consider for definiteness the case with ni = 77-2 = 0, 77-3 7^ 0. If ria+a = the space 
configuration has only two nests, iy^ea, see (6.88a). The structure is the same as in 
Figure 3: there are two nests. The equal probability the particle to be in one of them is a 
consequence of proposition 6.4 

The other opportunity is na+3 = 1, i.e. the state is = \p; 0,0,77-3, la+3) ••)• The 
space configuration is similar to the one on Figure 5. In cylindrical coordinates (p, a, z) 
the nests consist of all points with 

p = V2, < a < 27r, z = ±^^3 + 1. (7.52) 
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It is straightforward to verify that ^'3^' = and therefore 

gic^Ss^ = (7.53) 

i.e, the state is invariant under any rotations around 2;— axes. Let P(-\/2, iV'^s + 1) 
be the probabihty density the particle to be in the nest (-\/2, iV'^s + !)• Then from the 
rotation invariance around z one concludes that the probability density does not depend on 
whereas the parity invariance yields in addition that P(-\/2, ±-\/n3~+T) is independent 
on the sign in front of z. Then from the normalization condition for the probability density 
one finds P = l/\/327r. 

Finally we consider the Class I configurations. Then any point from a sphere with a 
radius y^3?v+3 is a nest for the ctthe particle. 

Proposition 7.5. Let the system he in a state \I/ = \p; 0, 0, 0; 724, ?1q,_(_3, ..) from Class I. 
Then the ath particle will occupy with equal probability any one of the allowed nests. 

Proof. If ria+s = the ctth particle is "sitting" on the center of mass with probability 1. 
So assume that na+3 = 1. It is straightforward to verify that Sf.'^f = for any k = 1, 2, 3. 
Therefore 

6^'^^'=* = *, A;= 1,2,3, (7.54) 

i.e, the state \E' is invariant under rotations around x—, y— or 2— axes. Hence it is invariant 
under an arbitrary rotation. Since any point a from this sphere can be moved by an 
appropriate rotation onto any other point b also from this sphere, the probability density 
is one and the same for any point on the sphere with radius -\/3- Then the normalization 
condition yields -P(a) = l/127r for any nest a of the sphere. 

Appendix A: Proof of Proposition 6.5. 

The transformation relations of the operators {raij'f^a2,'f^a3) under global rotations g are 
the same as for the frame vectors: 

3 

r{g)ai^U{g)raiU{g)~^ ^^Tajgji, i = 1,2,3, (A.l) 
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Proof. Since the results to be proved are the same for any of the particles, in what follows 
we skip the subscript a. We consider first a rotation ^(63, <^), namely a rotation around 
63 on angle (p. Then 

r{g{es, cp))^ = V{g{es, cp))rj = e-'^^H^e'^^^ , j = 1, 2, 3. (^.2) 

Since [<S'3,f3] = 0, the above relation yields 

f{9{e3,(fi))3 = f3- 

The computation of f{g{e3j(^))j for j = 1,2, is not that simple. Introduce first the 
eigenvectors of S3 (the weight vectors of the Cartan subalgebra), namely 

r+ = fi + ir2, and f_ = fi — zf2. (^-4) 

Then 

[S3,r+] = r+, and [83,7-^] = -r-. (A5) 

Next we compute e~'^'^^^f±e^'^^^ using the Backer-Campbell-Hausdorf formula (as given in 
[56] , which in this case reads: 

where the multiple commutator is defined as 

[A,B]^^^ = [A,B]=AB-BA, (A7) 
[AS](o) = S. 

Then 

H^^3,r±](fc) = (W)'=f±. {A.8) 
Inserting (A.8) in (A.6) one obtains: 

00 



k=0 



00 ( ■ \k 



(A.9) 



fc=0 
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And finally a replacement in (A. 9) of f± with fi, f2 according to (A. 4) yields: 

r{g{e^, (p))i = e~^'^^^rie'^'^^^ = fi cos (p + r2 sin (f, 

f{g{e3, (fi))2 = e~"^'^3f2e*'^'^3 = -fi sin (fi + f2 cos (fi, (A. 10) 

r{g{es,(fi))s = e-'^^'he"'^'=h 
The last result can be written in a compact form, 

3 

rigies, 'fi))j = V{g{e3, (fi))rj = hgies, (fi)k,j = {rg{e3, (fi))j, (^.11) 

k=l 

In the derivation of (A. 11) we have used only the commutation relations [S^, fj] which 
are invariant under cyclic change 

1^3, 2^1, 3^2. (A.12) 
Therefore Eqs. (A. 10) remain true under the change (A. 12): 

f(5r(e2, <^))3 = e-*'^^2 ^36**^^2 = cos (fi + h sin (fi, 
r{g{e2, (fi))i — e~'^'^^^rie'^'^^^ = — fa sin + fi cos 
f(^(e2,^))2 = e-^'^^V2e^'"^^ =f2, 

which yields 



3 

r{g{e2, (fi))j = "^^(^(62, (fi))rj = ^ h9{e2, (fi)k,j = {rg{e2, (fi))j. (A.13) 

/c=l 

Similarly, from (A. 12) and (A. 13) we derive 

r{g{Gu^))j = V{g{ei,(fi))rj = ^rkg{ei,(fi)k,j = {rg{ei,(fi))j, (^.14) 

fe=i 

Equations (A. 11), (A. 13), (A. 14) can be unified: 

3 

r{g{ei, (fii))j = V{g{ei, (fii))rj = ^ hgi^i, (fii)k,j = {rg{ei, (fi))j. {A.15) 
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We recall that V{g) gives a representation of SO (3). Therefore 
rig)i = y{g)h = V{g{e3,a)g{e2, P)gie3,^))ri 
= V{g(e3, a))V(gie2, PWigie^, 

= Yl ^idih, a)ikg{h, P)kjg{h, (A.16) 
= Yrig{e3,a)g{e2,/3)g{e3,'y))ii = ^figia, /3.'-f)ii = {rg)i. 



64 



The above result holds for any particle. In particular for the ath particle one has: 

3 

r{g)ai = r'»j9ji = {rag)i- (^-17) 

This completes the proof. 

Appendix B: Proof of Proposition 6.6 

// the system is in a basis state \p; n) from Class II, then the nests (6.83) are the only nests 
for the ath particle. 

Proof. We know from corollary 6.2 that if the dispersion Disp(f(5f)^ k)\p;n} vanishes for a 
certain g and for all A; = 1, 2, 3 then the set 

r{\p; n)), a, g) = {r{g)a,Mg)i + r{g)a,2e{g)2 + r{g)a,^e{g):i}. {B.l) 

determines admissible places, i.e., nests for the athe particle. We proceed to prove that 
the nests (B.l) coincide with those in (6.83) 

The first task to solve is to determine all 3 x 3 orthogonal matrices g for which the 
dispersion D{f{g)'^ k)\p;n) vanishes. According to (6.60) we have to find all solutions of the 
equation 

glk92ki'^nin2 + ni + 712) + glkglkC^nms + m + ns) 

{B.2) 

+ 93k92ki'^'n'2n3 + n2 + ns) = 0, /c = l,2,3, a = l,...,iV. 
where the unknown are the matrix elements of g. Since for all states from Class II 

2n2n3 + 77.2 + ns ^ 0, 2nins + ni + ns ^ 0, 2nin2 + ni + n2 ^ 0, (-^-3) 

the problem reduces to determine all solutions of the equations 

9lk9lk = 0, 9W3k = 0, 9lk9lk = 0, A; =1,2,3. (S.4) 

Here are the main steps in solving the problem. 

Find first all matrices which satisfy the restriction ^235'33 = 0. It yields two classes of 
matrices: 

Class 1. All matrices ^(a,/?, 7) with a, 7 being arbitrary and /? = 0,7r/2,7r, 
Class 2. All matrices g{a, 7) with 7 being arbitrary and a = 0, tt. 
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The conditions ^21^31 = and ^22^32 = do not lead to additional restrictions on 
Class 1. So we have: 

IA. All matrices g{a, (3 = 0, 7), which is a rotation of angle a + 7 about z axes. 

IB. All matrices g{a, /3 = tt, 7), which is a rotation of angle a — 7 about z axes. 
Ic. All matrices 5f(a, /? = 7r/2, 7). 

The equations g\f^g1j^ = 0, /c = 1, 2, 3 put additional restrictions only on the class Ic: 

IC. All matrices g{a, (3 = 7r/2, 7) with 7 = 0, 7r/2, tt, 3p/2. 

For further use we collect part of the results obtained so far. 

Corollary B.l. All solutions of the equations g2k93k = and gitgsk = 0, k = 1,2,3, 
are given with the g{a,P,'y) matrices from the subclasses lA, IB, IC, defined above (see 
below: the Class 2 does not contain new solutions). 

Finally the equations g\j^g2k = 0) /c = 1, 2, 3 lead to the following solutions: 

lAa 

(cos(a + 7), — sin(Q; + 7), 
sin(Q: + 7), cos(q; + 7), |, a + 7 = 0, 7r/2, tt, 37r/2, (S.5a) 
0, 1 

IBa. 

( — cos(q; — 7), — sin(Q; — 7), 
— sin(Q; — 7), cos(q; — 7), | , a — 7 = 0, 7r/2, tt, 3p/2. (B.bb) 
0, -1 

ICa. 

( — sin a sin 7, —sin a cos 7, cos a 
cosQ;sin7, cosq:cos7, sino: | , a, 7 = 0, 7r/2, tt, 3p/2. 
— cos 7, sin 7 

{B.5c) 

The conditions g2kgzk = is satisfied by the following matrices from Class 2: 

Class 2A: all g{a = 0, tt, /3 = 0, 7r/2, tt, 7 = 0, 7r/2, tt, 37r/2), 

Class 2B: aU g{a = 0, tt, /? = 0, tt, 7). 

Clearly the solutions from the Class 2A and Class 2B are particular cases of solutions 
from the classes lA, IB and IC and therefore we do not consider them anymore. 
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Thus if the system is in a state |p;n), for which conditions (B.3) holds, then the 
classes lAa, IBa and ICa determine all g matrices for which the dispersion of f{g)1^^. along 
e{g)ki k = 1, 2, 3 vanishes. Therefore for any such g the set 

3 

T{\p; n),a,g) = Y^ r{g)a,Mg)k (5.6) 

k=l 

determines a set of nests for the ath particle, where T^g)"^^ is an eigenvalue of f{gY^^, on 

^(f)a,fc =^lfe(^l+^a+3) +^ife(^2+n«+3) +^3fc("3 + na+3), /c = 1, 2, 3. {B .7) 

and r(£f)a,fe = ±^r{g)l ,^. 

At this place we shall make use of the following property of the matrices (B.2): 

Corollary B.2 Let g be any matrix from the classes lAa,lBa or ICa. Then each row and 
each column of g consist of two zeros and one number ±1. 

Then in view of the above property only one term in the RHS of (B.7) survives, 

rigfcck = +na+3) =^ r{g)a,k = ±^/n~^fri^ = raj^, {B.S) 

where {ji, j2, Js} is a permutation of {1, 2, 3} Similarly 

= Gj^gj^,k, {B.9) 

where gji^,k = 1 or —1. Therefore 

3 

r(|p; n)), a, g) = { r{g)a,ke{g)k} = ^V^j+^c,+3 = r(|p; n), a), (S.IO) 
k j=i 

see (6.24). This completes the proof. 



Appendix C: Proof of Proposition 6.8. 

Below we prove the third part of proposition 6.8, namely we show that the nests of the 
ath particle, whenever the system is in the state |p; 0, 0, 77,3, .., la+3, ..) are (6.91). The rest 
namely the equations (6.89) and (6.90)are proved in a similar way. 
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Consider first the solution lA , see corollary B.l. Without loss of generality we set 
= 7 = 0, leaving a to be arbitrary, i.e., we have 

cos a, —sin a, 
g{a, j3 — Q,'^ — Q) = | sin a, cos a, |, {C.l) 

0, 1 

which is a rotation about the 2;— axes on angle a. Then (6.87)yields the following eigen- 
values r{g)l^f. of rigfaj, on |p; 0, 0, na, n^+s, ..): 

r{9)li = l. r(^)a2 = l, K^)a3 = ^3 + 1. {C .2) 

Equation (6.68)teUs us that the nests of the ctth particle corresponding to g{a,P = 0, 7 = 
0), see (C.l), whenever the system is in the state \p; 0, 0,77-3, •••) la+3) ••) are: 

T{\p; 0, 0, ns, l^+s, ..), g{a, /3 = 0, 7 = 0), a) = ^ie{g)i + ^2e{g)2 + CsV^s + le(fi')3, 

(as) 

where ^1, ^2, 6 = ±1- 

In the initial basis e = ge{g) the nests (C.3) read 

r(|p;0,0,n3,...,lQ+3,..), 51(0;, /3 = 0,7 = 0), a) = {(^1 cosa - ^2sina)ei 

(C.4) 

+ (^1 sin a + ^2 cos a)e2 + C3\/«'3 + 1^3 | a e R}. 

The first impression might be that the nests corresponding to different choices of 
Ci) C25 Cs = are different. This is however not the case. With elementary considerations 
one shows that the different choices of ^1, ^2 lead to one and the same collection of nests. 
If for instance ^1 = 1, ^2 = 1, then the RHS of (C.4) reads: 

(cos a — sin a)ei + (sin a + cos a)e2 + ^sV^s + le(5r)3. 

Replacing in {C.5)a with a + 7r/2 one relabels the nests but does not change the collection 
of them. After this substitution (C.5)reads 

(— cos a — sin a)ei + (— sin a + cos 0)62 + ^sVns + le(^)3. (C'.6) 

which corresponds to the choice ^1 = —1, ^2 = 1 in (C4). Hence the choices ^1 = 1, ^2 = 1 
and ^1 = —1, ^2 = 1 describe one and the same collection of nests. In a similar way one 
concludes that also the other two choices ^1 = 1, ^2 = —1 and ^1 = —1, ^2 = —1 give the 
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same nests as the choice = 1, ^2 = 1- Thus the nests of the ath particle corresponding 
to solution lA, namely to the matrix (C.l)read: 



r(|p;0,0,n3, ...,la+3, ..), 5f(Q;,/3 = 0,7 = 0), a) = {(cosa - sinQ;)ei 
+ (sincK + cos a)e2 + ^sV'^s + lea | a e R, ^3 = ±1, }. 



iC.7) 



It takes some time to show that the nests corresponding to the IB and IC solutions 
of (6.85) whenever the system is in the state |p; 0, 0, 713, 1q,+3, ..) coincide with (C.7). 
Therefore the nests of the ctth particle corresponding to the state |p; 0, 0, n3, 1q_|_3, ..) 
are 

r(|p;0,0,n3, ...,la+3, ••), a) = {(cosa - sinQ!)ei 

(C.8) 

+ (sin a + cosQ;)e2 + ^363 | a G R, ^3 = ±1, }. 
This proves the third part (6.91 )of proposition 6.8. 

For an illustration set ris = 3. The corresponding nests (C.8) of ath particle in the 
state |p; 0, 0, 83, 1q,+3, ..) are indicated symbolically on Figure 3 as small black balls. 
There are infinitely many of them situated on two circles with radius y/2 around the 
2;— axes, which are on a distance 2 above and below the xOy plane. 



Fi^gure 5 




The space distribution of the nests of the ath particle, corresponding to an arbitrary 
state |p; 0, 0, 77-3, 1q;+3, ••), is similar: there are infinitely many nests situated on two 
circles with radius -\/2 around the 2;— axes, which are on a distance y/n^'+T above and 
below the xOy plane. 
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In a similar way one derives that the nests of the rest of the basis states from Class 
III are (6.89) and (6.90). 

Appendix D: Angular momentum structure of V{N,p) for p=l,2,3. 

The case p=l. According to conclusion (7.21) the state space V{N,p = 1) can be 
represented as follows: 

V{N,p = 1) = ViiN,p = 1, nb = 0) ® V2{N,p = 1, n/ = 1), 
®Vi{N,p=l,nb = l)^V2{N,p = l,nf = 0). 

In this case Vi{N,p = l,nf, = 0) and Vi{N,p = l,n(, = 1) are irreducible with respect to 
the rotation group. Therefore 

V{N,p = 1) = Vi{N,p = l^nt = O,S = O)0 V2{N,p - 1, n/ - 1), 

e Vi{N,p = l,n6 = 1, 5 = 1) «) V2{N,p= l,n/ = 0). {D.l) 

The case p=2. Again from (7.21) 

V{N,p = 2) = e(iV - 2)Vi{N,p^ 2, rife = 0) ® V2{N,p = 2, n/ = 2) 
© Vi{N,p = 2, nfo = 1) ® V2{N,p = 2,n/ = 1) 
© Vi{N,p = 2,n6 = 2) (8) V2{N,p = 2,n/ = 0.) 

This time the last subspace above is S'0(3) reducible. Therefore for the angular momentum 
content of V{N,p = 2) we obtain: 

V{N,p = 2) = e(iV - 2)Vi{N,p = 2,n;, = 0, 5 = 0) ® V2{N,p = 2, n/ = 2) 
©V^i(7V,p = 2,n6 = 1,5= l)(g)V^2(A^,P = 2,n/ = 1) 
© Vi{N, p = 2,nb = 2,S = 2)®V2{N,p = 2, n/ = 0) 
©Vi(iV,p = 2,n6 = 2,5 = 0)®y2(iV,P = 2,n/ = 0). (D.2) 

The case p=3. This time 
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V{N,p = 3) = Q{N - 3)Vi(7V,p = 3, = ^) ®V2{N,p = ^,nf = 3) 
e Q{N - 2)Vi{N,p = 3, 716 = 1) (8) V2iN,p= 3, n/ = 2) 

®Vi{N,p = 3,nb = 2)(g)V2iN,p = 3,nf = 1) 
® Vi(iV,p = 3, nfe = 3) ® V2(iV,p = 3,n/ = 0). 

The decomposition of each Vi into irreducible SO (3) modules yields (see (7.24): 

V(N,p = 3) = e(N - 3)Vi(N,p = 3, nfe = 0, 5 = 0) (8) V2{N,p = 3, n/ = 3) (D.3a) 
© e(iV - 2)Vi{N,p = 3, n6 = 1, >5 = 1) ® V^2(iV,P = 3,n/ = 2) (D.36) 
© Vi(iV,p = 3,n6 = 2, 5 = 2) O "l/2(iV,P = 3, n/ = 1) (i:».3c) 
© Vi(N,p = 3, n6 = 2, = 0) F2(iV,p = 3, n/ = 1) (r>.3d) 
© ViiN,p = 3,nb = 3,S = 3)^V2iN,p = 3,nf = 0) (L».3e) 
©Vi(iV,p = 3,nb = 3,5= l)®F2(iV,P = 3,n/ = 0) {D.3f) 

Since AT is the same for all basis vectors, instead of 1 1 AT, p, 77,5, /S, ^s, 77-4, njv+s)) we 
write I \p, 7i5, S, S3,n4, ...,nN^s))). The angular momentum basis vectors in each subspace 
(D.3a) - (D.3f), expressed via the reduced basis and via the initial basis, read 



D.3a) : 



{D.3b) : 



{D.3c) : 



\p, nb = 0,S = 0,S3 = 0; n^, ..)) = \p; 0, 0, 0; 714, ••) = \p; 0, 0, 0; n^, ■■), N ^ 1, 2; 



\p, Ub = 1,S = 1,83 = 1; 714, ••)) = -\p; 0, 1, 0; 714, ••) 
= --^ b; 1,0,0; 714, ••> - -^|p;0,l,0;7i4,..), 

\p, 716 = 1, 5 = 1, 53 = 0; ?i4, ..)) = \p; 0, 0, 1; 714, ..) = \p; 0, 0, 1; 714, ..); 

\p, Ub^l^S = 1,83 = -1; 714, ..)) = \p; 1, 0, 0; 774, ..) 

1 i 
= ;^b;l,0,0;7i4, ..) - --^|p;0,l,0;7i4, ..), N 1 



\p, 7ib = 2, 5 = 2, 53 = 2; 714, ••)) = \p; 0, 2, 0; 714, ••) 

1 z 1 

= 2b;2,0,0;7i4, ..) + --^b;l,l,0;7i4, ..) - 0, 2, 0; 714, ••); 

\p, rib = 2, S = 2,83 = 1; n4, ..)) = -\p; 0, 1, 1; 714, ..) 
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1 \ ^' I 

-;y||p;l,0, l;n4, ..) - -^Ip; 0, 1, 1; n4, ..); 



|p, n6 = 2, 5 = 2, ^3 = 0; n4, ..)) 



2 1 
- |p; 0, 0, 2; 714, ..) - -^b; 1, 1, 0; 714, ..) 



2 11 
-|p;0,0,2;n4, ..) - -^b; 2, 0, 0; 714, ..) - -^b; 0, 2, 0; 714, ..); 

I |p, rife = 2, 5 = 2, ^3 = -1; 714, ..)) = |p; 1, 0, 1; 714, ..) 

1 \ ^ I 

= ;^b;l,0, l;7i4, ..) - --^|p;0,l,l;7i4, ..),; 

I \p, 7ib = 2, 5 = 2, 53 = -2; 714, ..)) = \p; 2, 0, 0; 714, ..) 

1 i 1 

= 2b;2,0,0;7i4, ..) - ---^b;l,l,0;7i4, ..) + -jp; 0, 2, 0; 714, ••); 

(D.3d) : I |p, 716 = 2, 5 = 0, ^3 = 0; 714, ••)) = ^ Ip'^ 0' 0' 2; n^, ..) + ^\p-l,l, 0; 714, ••) 
= -^b; 0, 0, 2; 714, ••) + ^ b; 2, 0, 0; 714, ••) + ^ \p; 0, 2, 0; 714, ••>; 



{D.3e) : |b,7i6 = 3, 5 = 3,53 = 3; 714, ..)) = -\p; 0,3,0; n^, ..) = -^b; 3, 0, 0; 714, ••) 

— ^ b; 2, 1, 0, 7i4, ..) + — b; i> 2, 0, 714, ..) — —\p; 0, 3, 0, 714, ..); 

I \p, 3,S = 3,Ss = 2; n4, ..)) = b; 0, 2, 1; 714, ..) 

1 i I 

= 2b;2,0,l;7i4, ..) + -^\p;l,l,l;n4,..) - -b; 0, 2, 1; 714, ..); 

1 [4 
\\p,nb = 3,5 = 3,^3 = l;7i4, ..)) = -^b; 1> 2, 0; 714, ••) - y -b; 0, 1, 2; 714, ••) 

= b; 3, 0, 0; 7i4, ..) + b; 2,1,0, 714, ..) + b; 2, o; 714, ..) 

iv^i V V2, V v^i 

- ^^b;o>3,o,7i4, ..) - --^b;i>o>2;7i4, ..) + i;^b;0>i>2;7i4, ..); 

\\p,nb^3,S = 3,^3 = 0;7i4,..)) = -^^b; 0^ 0^ 3; 714, ••) - ^^b; 1, 1, 1; "4, ••) 

2" / 3 

-b;0,0,3;7i4, ..) - ]l — b;0:2, l;7i4, ••) - d — b; 2, 0, 1; 714, ..); 

1 [a 
\\p,nb = 3,5 = 3,^3 = -l;7i4, ••)) = ---^b;2, 1,0; 714, ..) + y -b; l,0,2;n4, 

y/g iA/2 \/2 

= -^;y|b;3,0,0;7i4, ..) + ^-^b; 2, 1, 0, 714, ••) - ^--^b;l>2,0;7i4, ..) 
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+ ^^b;0,3,0,n4, ..) + --^|p;l,0,2;n4, ..) - i-^jp; 0, 1, 2; n4, ..); 

\p, rib = ^,S = 3,83 = -2; n4, ..)) = |p; 2, 0, 1; 714, ..) 

1 i 1 

= 2b;2,0,l;n4, ..) - 1, 1, 1; 714, ..) + 0, 2, 1; 714, ..); 

V2 

|p, ri5 = 3, 5 = 3, ^3 = -3; 714, ..)) = |p; 3, 0, 0; 77,4, ••) = \p; 3, 0, 0; 714, ••) 
i\/6, , \/6, , 1^2, 

— ^|p;2, 1,0,714, ..) - 1,2, 0,714, ••) + 0,3,0,714, ••); 



{D.3f) : I |p, 716 = 3, 5 = 1, 53 = 1; 714, ..)) = b; 0, 1, 2; ^4, ••) - ^ ^ b; 1,2,0; 7^4, ••) 

3 i 1 

— |p; 3, 0, 0; n^, ..) 2, 1, 0, 774, ••) i=\p; 1, 2, 0; 774, ••) 

-LU vlO vlO 



/ 3 1 i 

+ TT^b^O' 3:0:^^4, ••) 7=|p; l,0,2;7i4, ..) + ^=|p; 0, 1, 2; 714, ••); 

V -LU vlO vlO 



\p,nb = 3,S = 1,83 = 0;n4,..)) = y 0, 0, 3; 714, ••) + y ^b; 1^ 1^ 1; ^4, ••) 

/3 1 1 

0, 0, 3; 714, ..) + -i=\p; 2, 0, 1; 714, ••) + -i=\p; 0, 2, 1; 714, ••); 



1 /4 
\p,nb = 3,S =1,83 = -l;7i4, ..)) = ;^b;l>0>2;7i4, ••) + y 2, 1, 0; 714, ..) 

1 i / 3 

, ^ b; 1> 0> 2; 7i4, ..) 7= |p; 0, 1, 2; 714, ••) + A/ ttt b; 3, 0, 0; 714, ..) 

ID vlO V iU 



i 1 / 3 

, ^ b; 2, 1, 0; 7i4, ..) + ^= \p; 1, 2, 0; 714, ••) - i\ — \p; 0, 3, 0, 714, ..)• 
lU vlO V -LU 



i /4 1 \/3 
2' V ^' V2' 2 
o^/b 1 

cv^' \/lO' V 10 



Appendix E: Eigenstates and eigenvalues of the momentum operators 

The A;th projection of the momentum operator for ath particle read (see (6.9b)) 



Pak{t) = ie(^Ek,a+3e''^' - Ea+3,ke-'"^'). {E.l) 
In the next proposition we write down the eigenstates of Pak{t) and their eigenvalues. 
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Proposition E.l. The eigenvectors of the momentum operator p a, k-, k = 1,2,3, read: 

a. Eigenvalue : 

«;°fe(..,0fe,..,0a+3,..) = b;..,Ofc, ...,0a+3,-), {E.2) 

b. Eigenvalues ±-s/nfe (^fc 7^ 0) : 

± i£(_l)ni + ...+n„+.g-iea.t|^. ^ _ ^ ^^^^^ > (^ 3) 

T/ie inverse to (E.3) relations take the form (10370): 

\p\ .., nfe, .., Oa+3, ••) = («^afc(--> Oa+3, ••) + ^^afcl"' -, Oa+3, ••)) , 

|p; .., nfc - 1, .., la+3, ..) = -l=e(-l)(-i+-+"«+^)e--* {EA) 

The same equations in a compact form: 

|p;..,nfc,,..,n«+3,..) = ^(_i)(-i+-+-«+^+i)-«+3e-"«+3c.t 

(^afc(--' ^fc + ^a+3> Oa+3, ••) + (-l)''"+'«^afc(--' + ^ia+3, -, Oa+3, ••)) , i^-^) 

Note that the eigenvectors of the momentum operators ^ and the eigenvectors of 
the position operators fa^ki corresponding to zeroth eigenvalues coincide, 

Wafc(", Ofc, Oa+3,..) = •^afcl-, Ofe, Oa+3,..) = |p; ••, Ofc, Oa+3, ••)• (^-6) 

Therefore fa,fe and Pa,fc commute on the subspace spanned by aU states \p] .., 0^, .., Oq+s, ..) 



8. Concluding remarks 

We have studied the properties of A?"— particle noncanonical harmonic oscillator, consider- 
ing it as a Wigner quantum system with the additional requirement the position and the 
momentum operators of the oscillating particles to be odd operators, generating the Lie 
superalgebra sl{Z\N). 
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The idea for such a requirement is a natural one if one takes into account that the 
canonical PM-operators generate also a representation of a LS, but from the class B, namely 
the orthosymplectic LS osp{l\QN). We should admit however that such an assumption is of 
pure mathematical origin and as such the sZ(3|A^) oscillator is essentially a mathematical 
model. Nevertheless it is surprising to see how rich is the idea of Wigner to relax the 
postulates of QM replacing the postulate about the CCRs with the requirement both the 
Heisenberg and the Hamiltonian's equations to be fulfilled simultaneously. 

We see that despite of the circumstance that the equations of motion (1.5), the Heisen- 
berg equations (1.6) and the Hamiltonian (1.1) are formally the same as for N free oscil- 
lators, the properties of the sl{3\N) oscillator are very different from those of the corre- 
sponding canonical such oscillator. 

On the first place Conclusion 5.9 tell us that there exist strong space correlations 
between the particles. In particular all particles from Class II, which have one and the 
same fermionic coordinates "share" 8 common nests independently on the number N of the 
"inhabitants", the oscillating particles. Clearly these correlations are of statistical origin. 

Secondly, and this is an essentially new result, there exists even stronger statistical 
correlation between the angular momenta of the particles: the components of the angular 
momentum of all N particles Mai, -Mq,2, M^s coincide, they do not depend on the label of 
the particle a, see (3.30). Consequently all particles have one and the same angular momen- 
tum. For instance if the system is in a reduce basis state | \N,p, b, S, Ss, n^, n^, njv+a)), 
see (7.28), then all particles have one and the same angular momentum S and one and the 
same projection along 2;— axes S^. 

Another property to mention is the space structure of the basis states \p; n) . Typically 
each such state corresponds to a picture when each oscillating particle is measured to 
occupy with equal probability only finite number of points, as a rule the eight vertices of 
a parallelepiped. As a result the entire oscillator is confined in the coordinate (and the 
momentum) space, it is "locked" within a sphere with a finite radius. This property is in 
the origin of the relations Af < y^, Ap < y^. Therefore in the limit ^ — > the entire 
system collapses into a point. 

We should not forget to point out again that all our conclusions and results hold only 
for the Hamiltonian (1.1). Is the idea of Wigner applicable for another Hamiltonians and 
what are their predictions is an open question. A first example in this directions was given 
in [20] for a magnetic dipole precessing in magnetic field 
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The results obtained in the present paper are based on a pure class of irreducible rep- 
resentations of sl{3\N), which are numbered by only one positive integer p. The reason for 
such a choice stems from the observation that explicit expressions for all finite-dimensional 
representations of the LS s/(3|A^) does not exists so far. 

Therefore a natural next step would be to extend the class of representations of the 
underlying superalgebra. The simplest way to do this is to consider the realizations of the 
superalgebra s/(3| A^) in the Fock space of three pairs of Fermi CAOs (even generators) and 
N pairs of Bose CAOs (odd generators), where again the Bose CAOs anticommute with 
the Fermi operators. This realization will put strong limitations on the angular momentum 
of the entire system. Another possibility is to consider the Holstein-Primakoff realization 
of gl{3\N + 1) [N 119]. In this case the Bose operators are even operators and the Fermi 
CAOs are odd and Bose operators commute with Fermi CAOs. In such a case the Fock 
space wouldn't be anymore an irreducible module of an orthosymplectic LS but of the 
more familiar gl{3\N + 1). 

As already mentioned, the main problem arising in the context of WQS is to determine 
the common solutions of the equations of motion (1.5) and of the Heisenberg equations 
(1.6), which satisfy the defining postulates (PI) - (P6). Stated in this way, the problem 
does not require the PM-operators to be elements of a Lie superalgebra, of a Lie algebra, 
or of any another algebraic structure. From this point of view the LS fif/(3|A^ + 1) was only 
a tool to find at least some solutions of the problem. And these solutions turned to predict 
WQS with interesting properties. 

Are there any indications that WQSs exist in nature? Can they be of real interest in 
physics? In this relation we mention that recently the finite-level quantum systems become 
of great interest in quantum computing. And such are the WQSs. The properties of the 
WQSs reported in the present paper resemble also the artificial atoms arising in condensed 
matter physics [61] and more generally various kinds of clusters (see, for instance [62]). Is 
there any deeper connection behind just a resemblance? That is what we would like to 
know too. 
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